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ABSTRACT 


A  computer  program  for  determining  the  natural  frequency  of  a  ship  hull  elastically 
connected  to  other  elastic  systems  and  to  sprung  masses  and  containing  concentrated  and 
hysteresis  damping  is  developed.  This  program  is  applied  to  the  mass-elastic  system  that 
represents  the  SSB(N)  598  GEORGE  WASHINGTON  in  detail.  The  influence  upon  the 
hull  response  to  propeller  and  hull  excitation  in  the  stern  of  variables  such  as  the  number 
of  mass  points  used  to  represent  the  hull,  the  omission  of  the  propulsion  sub-system,  the 
stiffness  between  the  propulsion  sub-system  and  the  primary  hull,  the  inclusion  of  the 
moment  of  inertia  of  the  hull  cross-sections,  the  treatment  of  water  inertias,  the  treatment 
of  sprung  masses  and  the  amount  of  hysteresis  damping  in  the  hull  is  shown. 


viii 
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SUMMARY 


STATEMENT  OF  THE  PROBLEM 

2 

A  previous  prediction  of  the  amplitude  of  longitudinal  and  vertical  bending 
vibration  excited  by  the  propeller  for  the  submarine  GEORGE  WASHINGTON,  SSB(N)  598, 
gave  results  that  agreed  moderately  well  with  tests  in  longitudinal  vibration  and  did  not 
agree  for  vertical  bending  vibration.  This  study  has  been  undertaken  to  explore  the  factors 
that  influence  the  structural  response  of  the  hull  and  eventually  (in  a  subsequent  report)  to 
see  whether  it  is  possible  to  make  a  more  accurate  prediction  of  the  amplitudes  of  propeller 
excited  hull  vibration 


WORK  SUMMARY 

1.  A  computing  machine  program  that  is  sufficiently  flexible  to  deal  with  parallel 
sub-systems  (e.  g ,  the  propulsion  system),  sprung  masses  and  concentrated  and 
hysteresis  damping  in  axial  and  bending  vibration  has  been  developed. 

2.  The  mass-elastic  characteristics  of  the  SSB(N)  598  GEORGE  WASHINGTON 
have  been  determined  in  considerable  detail. 

3.  The  effects  of  various  approximations  such  as  omitting  the  propulsion  sub-system, 
treating  sprung  masses  as  rigidly  connected,  neglecting  rotary  inertia,  repre¬ 
senting  the  hull  with,  a  limited  number  of  mass  points,  spacing  the  mass  points 
uniformly  rather  than  to  conform  to  hull  mass  and  stiffness  variations,  have  been 
investigated 

4.  The  effects  of  various  amounts  of  hysteresis  damping  in  the  hull  have  been 
investigated. 


ix 
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CONCLUSIONS 

1.  In  predicting  the  vibration  response  of  a  ship  in  axial  and  bending  modes,  it 

is  necessary  to  include  the  propulsion  system  (propellers  and  shafting)  as  an  elastic  sub-system 
in  the  calculations  if  reliable  results  are  to  be  obtained. 

2.  Hull  damping  not  only  reduces  the  vibration  amplitude  in  the  vicinity  of  the 
excitation  but  also  attenuates  the  relative  amplitude  at  locations  remote  from  the  point  of 
excitation. 

3.  In  order  to  obtain  reliable  predictions  of  vibration  response  of  a  ship  it  is 
necessary  to  define  the  hull  mass  and  stiffnesses  accurately.  A  less  accurate  definition  is 
required  if  only  hull  frequencies  are  desired.  Because  of  the  attenuation  caused  by  damping, 
the  accurate  definition  of  the  hull  is  particularly  important  in  the  regions  near  the  excitation. 

4.  On  the  submarines  the  effect  of  including  rotary  inertia  of  the  hull  cross 
sections  upon  the  response  of  the  ship  are  so  small  as  to  be  negligible. 

5.  When  a  careful  effort  has  been  expended  to  define  the  hull  mass  and  stiffnesses 
accurately,  it  is  not  advisable  to  use  only  a  few  mass  points  to  represent  the  hull.  The  com¬ 
putation  program  can  accept  many  points  at  a  negligible  extra  cost.  The  greater  number 

of  points  will  result  in  better  accuracy  and  a  better  definition  of  the  hull  motions  will  be 
obtained. 

6.  The  hull  response  is  sensitive  to  the  stiffness  of  the  connections  between  the 
propulsion  system  and  the  hull,  and  therefore  careful  attention  should  be  paid  to  these 
stiffnesses.  By  controlling  these  stiffnesses  to  locate  the  resonant  frequencies  between  the 
propulsion  sub-system  and  the  hull,  it  is  possible  to  influence  the  hull  vibration  amplitudes 
at  low  speeds  as  compared  with  those  at  high  speeds. 

7.  Although  the  total  weight  of  the  sprung  masses  is  small  compared  to  the  weight 
of  the  hull,  including  them  as  sprung  weights  rather  than  as  solidly  mounted  weights  has  an 
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important  effect  upon  the  hull  response. 

RECOMMENDATIONS 

If  reliable  estimates  of  the  response  of  a  ship  hull  to  vibratory  excitation  are  to 
be  obtained  it  is  necessary  to  determine  the  mass  and  elastic  characteristics  of  the  hull 
in  considerable  detail.  Particular  attention  should  be  paid  to  the  propulsion  sub-system 
and  its  connection  to  the  hull.  The  local  springing  of  hull  weights  must  be  included. 


xi 
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I  INTRODUCTION 

In  an  early  report^  CONESCO  outlined  a  program  for  predicting  the  amplitude  of 

propeller  excited  hull  vibration  and  presented  ways  of  estimating  the  excitation  arising 

2 

from  the  propellers  working  in  irregular  wakes.  A  subsequent  report  applied  this  program 

to  a  prediction  of  the  longitudinal  and  vertical  bending  propeller  excited  vibrations  on  the 

USS  GEORGE  WASHINGTON,  SSB(N)  For  this  purpose  the  excitation  was  computed 

by  the  methods  of  the  earlier  report  but  the  structural  response  of  the  hull  was  determined 

3 

from  existing  calculations  .  |t  was  considered  that  if  the  ship  resonant  frequencies  as 
measured  on  trials  and  the  amplitudes  of  vibration  on  the  flanks  of  the  resonances  agreed 
with  those  that  were  predicted,  that  the  method  of  prediction  would  be  validated.  The 
amplitude  of  vibration  at  the  resonances  would  depend  upon  more  detailed  considerations 
of  damping  in  the  calculations  and  more  experimental  comparison  between  the  damped 
calculations  and  trial  results  than  could  be  made  at  the  time. 

When  this  comparison  was  made  it  was  found  that  there  was  a  moderate  support  for 
the  calculation  method  as  applied  of  longitudinal  vibration  of  the  hull  since  the  predicted 
amplitudes  of  vibration  on  the  flanks  agree  moderately  well  with  measured  values  and  the 
calculated  and  measured  natural  frequencies  agreed  within  about  10%.  For  the  bending 
calculation  there  appeared  to  be  no  agreement.  However,  the  representation  of  the  hull 
in  the  axial  response  calculations  took  account  of  the  elasticity  of  the  connection 
between  the  propeller  and  the  hull  where  the  bending  calculation  did  not. 

In  an  effort  to  improve  the  accuracy  of  the  computation  process  and  to  study  the 
effects  of  damping  on  the  hull  vibration  amplitude,  a  thorough  study  of  the  structural  response 
of  the  hull  has  been  undertaken.  To  accomplish  *his  it  has  been  necessary  to  develop  a 
new  program  for  computing  hull  response  in  bending  which  is  sufficiently  adaptable  to 
include  not  only  the  mass  and  elasticity  of  the  elements  in  the  propulsion  system,  but 
also  the  elastically  supported  elements  throughout  the  ship.  A  procedure  for  computing 
the  response  of  a  ship  that  is  represented  in  the  most  general  way  of  a  beam  having 
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distributed  mass  and  stiffness  was  also  considered  but  because  of  complexity  in  the  pro¬ 
gramming  could  not  be  completed  and  is  therefore  not  reported  at  this  time. 

The  results  of  these  studies,  involving  (1)  a  more  precise  representation  of  the 
hull  mass  and  structure  (2)  the  inclusion  of  elastic  sub-systems  and  (3)  an  investigation 
of  the  effects  of  simplifications  and  assumptions  made  in  the  calculations,  are  presented 
in  this  unclassified  report.  Because  the  actual  vibration  of  the  ship  falls  in  the  realm  of 
classified  information,  the  comparison  of  calculated  and  experimental  amplitudes  of 
vibration  will  be  presented  in  a  subsequent  classified  report. 

II  METHOD  OF  SOLUTION 

The  method  of  computing  the  amplitudes  of  the  propeller  excited  hull  vibration 
falls  into  two  categories:  (1)  the  calculation  of  the  exciting  forces  and  (2)  the  calculation 
of  the  response  of  the  hull  to  these  exciting  forces.  This  is  essentially  the  standard  pro¬ 
cedure  followed  in  the  calculation  of  the  forced  vibrations  of  any  system  but  in  the  case 
of  ships  involves  many  difficulties.  The  flow  of  work  in  the  prediction  process  is  dia- 

2  t 

grammed  in  Figure  1.  The  prediction  of  the  forces  has  been  discussed  in  previous  reports.  ' 
This  report  is  concerned  solely  with  the  response  of  the  hull  to  these  forces. 

III  DEFINITION  OF  THE  MASS  AND  ELASTIC  PROPERTIES  OF  THE  SHIP 

|t  should  be  pointed  out  that  since  the  ship  under  consideration  is  a  submarine,  the 
definition  of  the  mass  and  stiffness  can  be  made  quite  exact.  The  requirement  of  adequate 
transverse  strength  results  in  heavy  structural  members  which  carry  their  loads  without 
buckling  or  local  resonances.  Furthermore,  the  shape  is  readily  adaptable  to  theoretical 
predictions  of  the  inertias  of  entrained  water.  It  should  not  be  assumed  however,  that 
even  the  submarine  can  be  defined  precisely.  The  flexibilities  and  water  inertias  associated 
with  the  control  surfaces  are  inadequately  handled  in  the  present  study  and  there  are 
questions  about  the  stiffness  of  the  superstructure  and  the  mass  of  the  free  flooding  water 
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SEQUENCE  OF  CALCULATIONS  FOR  ESTIMATING 
THE  AMPLITUDE  OF  PROPELLER  EXCITED  HULL  VIBRATION 


I kutnndv  Flow  Calculation)  I  I  riHull  Lined  1  .  I  MCnnfinuration! 


Figure  1  -  Calculation  Sequence 
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between  the  superstructure  and  the  hull. 

In  this  study  the  mass  distribution  and  the  stiffness  distribution  were  determined 
to  a  rather  high  degree  of  precision  and  then  plotted.  From  these  plotted  results  a  variety 
of  hull  representations  were  developed  which  seemed  to  fit  the  distributions  of  mass  and 
stiffness  well.  A  small  scale  plot  of  stiffnesses  is  given  in  Figure  2  and  of  weight  dis¬ 
tribution  in  Figure  3.  The  weight  moments  of  inertia  of  the  hull  are  presented  in  Figure  4. 
(Apparently  E.  B  division  feels  that  the  entrained  water  contributes  to  moment  of  inertia. 

We  feel  that  the  contribution,  as  to  the  axial  mass,  is  negligible).  The  weight  and 
stiffness  of  the  propulsion  system  in  longitudinal  vibration  is  given  in  Figure  5  and  in 
vertical  vibration  in  Figure  6.  The  characteristics  of  the  sprung  mounted  weights  and  their 
mountings  are  given  in  Table  I.  A  discussion  of  the  methods  by  which  these  quantities 
were  determined  and  the  factors  that  entered  into  their  determination  are  given  in  Appen¬ 
dices  A-H, 

,V  CALCULATION  PROCEDURE 

For  the  longitudinal  vibration  it  was  found  that  an  existing  Taylor  Model  Basin 
Code  was  adequate  for  handling  all  of  the  variables  that  it  was  desired  to  impose  upon 
the  problem  For  the  analysis  of  the  bending,  all  of  the  existing  codes  were  quite  inadequate. 
CONESCO  O'iginally  proposed  some  additions  by  which  the  existing  codes  might  be  made 
workable  but  the  addition  of  these  additional  complications  to  an  already  inefficient  code 
was  not  considered  desirable.  Dr.  Cuthill  (one  of  the  authors  of  this  report)  therefore 
prepared  an  entirely  new  code  which  is  adaptable  to  a  wide  variety  of  hull  representations. 
This  code  can  deal  not  only  with  the  bending  but  also  is  sufficiently  general  to  cover  the 
longitudinal  vibration  and  with  small  changes  can  deal  with  coupled  longitudinal  and 
vertical  bending  or  coupled  torsion  and  transverse  bending.  Because  it  was  developed  after 
the  hull  had  been  defined,  the  present  calculations  do  not  fully  exploit  the  potentials  of 
the  new  code  (specifically  in  the  treatment  of  control  surfaces).  A  presentation  of  the 
basis  for  this  code  and  the  way  it  is  used  is  given  in  Appendix  I.  This  general  bending 
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.Figure  2  -  SSB(N)  598  -  Axial,  Shear  and  Bending  Stiffness 
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Distance  from  Forward  Perpendicular  -  Feet 
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Figure  5-SSB(N)  598  Weightand  Stiffness  of  Propulsion  System  in  Longitudinal  Vibration 


Stiffness  of  Propulsion  System  in  Vertical  Vibration 
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Damping  is  located  between  mass  and  foundation  except  for  fainvater  and  stem  p.  tes  which  are  to  ground,  Fair- 
«-ater  and  stem  plane  damping  constant  is  proportional  to  ship  speed.  This  proportional  ity  is  represented  by  n,  the 
frequency  of  oscillation,  cps,  for  the  case  of  a  5-bladed  propeller. 


CONESCO,  INC.,  consultants  in  engineering  science 

response  code  carries  the  AML  problem  number  designation  840-277. 

V  DAMPING 

Damping  on  a  complex  structure  like  a  ship  is  difficult  to  define  —  in  fact,  it  can 
only  be  determined  by  experience.  However,  there  are  certain  types  of  damping  that  can 
be  calculated  by  theory  or  from  the  results  of  experiments  external  to  the  ship.  Among  these 
are  the  damping  of  the  propeller  in  longitudinal/orsional  and  transverse  motion,  the  aero¬ 
dynamic  damping  of  control  surfaces  and  the  damping  in  the  rubber  mounts  used  under  flexibly 
mounted  machinery.  Although  these  are  significant  sources  of  damping,  the  major  damping 
is  in  the  hull.  By  experiments  and  theory  it  is  known  that  this  damping  is  in  the  hull  itself 
not  in  the  water  and  is  thus  related  to  the  slope  of  the  deflection  curve,  not  the  at  solute 
amplitude,  i.e. ,  the  relative  amplitudes  between  adjacent  points.  In  the  damping  of 
structures  it  has  been  found  that  one  of  the  most  reasonable  ways  to  deal  with  this  type  of 
damping  is  to  tie  it  with  the  stress  in  the  structure.  |f  it  is  assumed  that  a  certain  percentage 
of  the  stress  energy  in  a  volume  of  stressed  material  is  dissipated  because  of  damping,  then 
it  is  possible  to  incorporate  the  effects  of  this  damping  in  the  calculations  by  the  simple 
expedient  of  making  the  modulus  of  elasticity  a  complex  term.  It  should  be  understood 
that  the  actual  damping  in  a  ship  is  much  larger  than  that  due  to  the  hysteresis  in  the 
material.  Most  of  the  damping  probably  comes  from  sliding  between  elements  of  structure 
and  cargo.  It  is  tied  to  stress  only  because  this  procedure  offers  a  good  basis  for  calculations, 
and  is  not  incompatible  with  experience.  The  actual  value  of  the  damping  constants  can 
only  be  determined  by  many  comparisons  between  carefully  calculated  amplitudes  of 
vibration  and  those  obtained  on  trials.  The  damping  constants  will  probably  vary  from  ship 
to  ship.  |t  will  be  noted  in  the  subsequent  results  that  the  damping  not  only  reduces  the 
amplitudes  ^  vibration  in  the  vicinity  of  the  exciting  force  but  also  reduces  the  transmission 
of  the  vibration  into  parts  of  the  structure  remote  from  the  excitation. 
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VI  REPRESENTATION  OF  THE  HULL  IN  THE  CALCULATIONS 

Because  the  hull  was  so  carefully  defined  and  because  a  computer  program  had 
been  developed  that  could  deal  with  a  wide  variety  of  hull  representations,  the  stage  was 
well  set  for  a  broad  exploration  of  many  factors  that  were  expected  to  influence  the  hull 
response. 

The  various  calculations  that  were  made  on  the  hull  of  the  SSB(N)  598  are  given  in 
Table  II.  Supporting  data  giving  the  constants  for  longitudinal  vibration  are  given  in 
Figures  6,  7  and  8.  For  the  bending  vibration  calculations,  the  arrangement  of  the  stations 
are  given  in  Figures  9-14  in  conjunction  with  Table  III.  The  location  of  the  stations  is 
given  in  Table  IV  The  specific  values  of  mass  and  stiffnesses  in  bending  vibration  are  given 
in  the  computer  input  sheets  in  Appendix  J. 

In  both  the  longitudinal  and  the  bending  vibrations  the  fi'St  computation  was  a 
repeat  of  the  earlier  one  made  in  Reference  3  (Cases  i  ard  12).  In  addition,  a  repeat  of 
the  calculations  with  a  10%  hysteresis  damping  was  made  in  orde'  to  compare  with  !ater 
calculations  (Cases  la  and  12a).  The  second  calculation  was  of  a  hull  represented  by  the 
same  number  and  location  of  stations  as  in  Reference  3,  but  with  each  station  assigned  the 
mass  and  stiffness  determined  by  the  more  detailed  anaiysis  of  the  hull.  (Cases  2  and  13), 

The  comparison  between  these  calculations  gives  an  indication  of  the  possible  increase  in 
accuracy  that  could  be  gained  by  the  laborious  defir>:»iOP  Qf  the  hull.  In  longitudinal 
vibration  an  exploration  was  made  of  *he  number  of  mosses  used  to  represen*  the  ship  (Cases  2, 
3  and  4);  of  *he  effects  of  hull  damping  (Cases  4  5,  6  and  7)  and  of  the  effects  of  changes 
in  the  difficult-*o-compute  stiffness  between  *he  p-opulsion  system  and  the  hull  (Cases  4,  8 
and  9).  The  final  two  runs  (Cases  10  and  i  ■)  we  e  mode  fo>  use  in  computing  the  hull 
vibration  and  indicate  the  difference  in  hull  'espcnse  for  fo'ces  on  *he  propeller  and  on  the 
hull. 

For  the  bending  vibrations  there  «ere  additional  foc*o's  to  be  investigated  and  so 
the  range  of  calculations  was  large’.  The  difference  in  hull  response  as  the  propeller  and 
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Text  continues  on  page  34 
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Vertical  "  Unit  vertical 

Bending  Figure  1 1  force  at 

Appendix  Propeller 
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OF  HULL  VIBRATION  CALCULATIONS  (continued) 
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TABLE  II  -  SCHEDULE  OF  HULL  VIBRATION  CALCULATION  (continued) 
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TABLE  II  -  SCHEDULE  OF  HULL  VIBRATION  CALCULATIONS  (continued) 
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TABLE  II  -  SCHEDULE  OF  HULL 
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The  sectional  water  inertia  values  are  reduced  to  account  for  bending  vibration  modes  by  the  value  J. 
See  Appendix  C. 
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Representation  Figure  12  -  Ship  Representation  for  Calculation 
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TABLE  III  -  35  MASS  REPRESENTATION  -  VARIATIONS  OF  VERTICAL  BENDING 

VIBRATION  CALCULATIONS 


Case  No. 

Water 

Inertia 

Hull  Damping 

Excitation 

Other 

17 

J  =  0.  87  for 
0<n  <3.  5 
=  0.  82  for 
3.5<n 

10%  Hysteresis 

Vertical  force 
at  Propeller 
(Sta.  56) 

18 

Standard 

J  determined 
for  sinusoidal 
mode  shape 

10%  Hysteresis 

Vertical  force 
at  Propeller 
(Sta.  56) 

19 

Standard 

No  hull  damping 

Standard 

20 

Standard 

4%  Hysteresis 

Standard 

21 

Standard 

25%  Hysteresis 

Standard 

22 

Standard 

Standard 

Standard 

Stiffened 
connections 
between  propul¬ 
sion  system  and 
hull 

23 

Standard 

Standard 

Standard 

Softened 
connections  be¬ 
tween  propulsion 
system  and  hull 

24 

Standard 

Standard 

Moment  about 
the  horizontal 
axis  at  the 
propeller 

25 

Standard 

Standard 

Vertical  force 
at  stern  (Sta,  54) 

26 

Standard 

Standard 

Vertical  force 
forward  of  the 
stern  (Sta.  51) 

30 

Standard 

Standard 

Standard 

Include  rotary 
inertia  of  the 
hull  cross-section 

Excitation  at  Propeller  (Sta.  39) 

10%  Hysteresis  Damping 

Figure  15  -  Ship  Representation  for  Calculation 
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(s)  indicates  sprung  mass 
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(s)  indicates  sprung  mass 
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its  connection  to  the  hull  is  brought  into  consideration  is  shown  by  a  comparison  between 
Cases  13  and  14  where  the  excitation  remains  on  the  hull  and  with  Case  15  where  the 
excitation  is  transferred  to  the  propeller.  The  influence  of  the  number  of  masses  used  to 
define  the  ship  is  shown  by  a  comparison  of  Cases  (5,  16  and  17.  A  comparison  between 
runs  17  and  18  indicates  the  effects  upon  the  hull  response  of  different  assumptions  for  the 
reduction  of  the  sectional  water  inertia  to  allow  for  the  modal  pattern  of  the  hull  vibration 
(See  Appendix  C).  The  effects  of  hull  damping  are  determined  through  Cases  18,  19,  20 
and  21.  The  effects  of  changes  in  the  stiffness  of  the  connections  between  the  propulsion 
system  and  the  hull  are  shown  by  Cases  18,  22  and  23.  Case  24  is  included  to  give  a  base 
for  estimating  the  hull  response  to  a  moment  at  the  propeller  and  Cases  25  and  26  to 
forces  acting  on  the  hull.  Case  27  is  included  for  comparison  with  Case  18  to  show  the 
influences  of  the  sprung  masses  (whose  gross  weight  is  very  small)  upon  the  hull  response 
and  Case  30  by  comparison  with  Case  18  shows  the  importance  of  including  the  rotary 
inertia  of  the  hull  cross  sections  in  making  calculations. 

VII  RESULTS  OF  CALCULATIONS 

The  results  of  the  calculations,  which  o'e  in  the  form  of  deflections  at  each  station 
for  a  range  of  frequencies  from  1  to  25  cps,  are  best  shown  by  curves  of  hull  motion  per 
pound  of  excitation  force  at  selected  positions  and  phase  angle  relative  to  the  excitation. 
These  curves  are  presented  in  Figures  16  through  60.  From  a  study  of  these  curves  it  is 
possible  to  determine  the  answers  to  some  of  the  questions  raised  in  the  previous  section 
of  this  report. 

1.  Longitudinal  Vibration: 

In  comparing  the  results  of  the  several  longitudinal  vibration  calculations  there  are 
several  lessons  to  be  noted  but  there  does  not  appear  to  be  any  trend  that  could  not  be  ex¬ 
pected.  Case  7(Figures  21  and  22)  is  the  standard  to  which  other  cases  are  compared. 

In  comparing  Case  1  with  Case  7,  the  first  noticeable  difference  is  the  larger 
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Frequency,  cpi 

Figure  16- Response  Curves  in  Longitudinal  Vibration 


Case  la.  Excitation  at  Propeller.  10%  Hysteresis  Damping 
Mass  and  Stiffness  taken  from  DTMB  Report  1464 
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Figure  17  -  Response  Curves  in  Longitudinal  Vibration 


Case  3.  Excitation  at  Propeller.  10%  Hysteresis  Damping 
24  Uniformly  Spaced  Masses,  Sprung  Masses,  Propulsion  System 


Figure  18  -  Retponte  Curvet  in  Longitudinal  Vibration 

I  Com  4.  Excitation  at  Propel  lor.  10%  Hystoretit  Damping 

I  16  Concentrated  Manet,  Sprung  Manet,  Propultion  System 
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Figure  19  -  Response  Curves  in  Longitudinal  Vibration 

Case  $.  Excitation  at  Propeller.  No  Hysteresis  Damping 
35  Concentrated  Mosses,  Sprung  Masses,  Propulsion  System 
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Figure  20  -  Response  Curves  in  Longitudinol  Vibration 

Cose  6.  Excitation  at  Propeller.  4%  Hysteresis  Damping 
35  Concentrated  Masses,  Sprung  Masses,  Propulsion  System 

I 
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Frequency,  cpt 


Figure  21  -  Response  Curves  in  Longitudinal  Vibration 

Case  7.  Excitation  at  Propeller.  10%  Hysteresis  Damping 
35  Concentrated  Masses,  Sprung  Mosses,  Propulsion  System 


Frequency,  cps 

‘Figur*  22  -  Longitudinal  Vibration  -  PhaM  Angle  Relative  to  Propeller  Excitation,  Cote  7 


Frequency,  cpt 


Figure  23  -  Response  Curvet  in  Longitudinal  Vibration 

Com  8.  Excitation  at  Propeller.  25%  Hyeiteretit  Damping 
35  Concentrated  Mauet,  Sprung  Mattes,  Propulsion  System 
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Frequency,  cps 


Figure  24  -  Retporae  Curvet  in  Longitudinal  Vibration 


Cate  9.  Excitation  at  Propeller.  10%  Hytteretl*  Damping 
35  Concentrated  Mattel,  Sprung  Matte*,  Propultion  Syttem,  Stiffnett  of  Connection 
between  Hull  and  Thrutt  Bearing  twice  that  of  Cate  7. 
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Frequency,  cps 
Figure  25- Response  Curves  in  Longitudinol  Vibration 
Case  10.  Excitation  at  Propeller.  10%  Hysteresis  Damping 


35  Concentrated  Masses,  Sprung  Masses,  Propulsion  System,  Stiffness  of  Connection 
between  Hull  and  Thrust  Bearing,  30%  of  that  in  Case  7. 
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Com  11.  Excitation  at  Stom.  10%  Hystorosis  Damping 
35  Concontratad  Mass**,  Sprung  Mono*,  Propulsion  System 
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amplitude  of  hull  motion  in  Case  1.  The  resonant  frequencies  occur  at  about  the  same 
frequencies  (except  that  the  propeller  in  Case  7  shows  a  weak  resonance  at  about  5  cycles 
that  is  not  present  in  Case  1.)  This  corresponds  to  a  mode  in  which  the  high  speed  parts 
of  the  system  (turbines  and  high  speed  gears)  vibrate  against  the  propeller  and  the  hull 
through  the  soft  longitudinal  stiffness  of  the  bull  gear  web.  The  differences  in  hull  amplitude 
can  be  explained  by  the  addition  of  propeller  damping  and  the  damping  of  the  sprung  masses 
to  the  general  hysteresis  damping  of  the  hull  in  Case  7. 

A  comparison  of  the  effect  of  the  number  of  stations  used  in  defining  the  hull  can 
be  obtained  by  considering  Cases  3  (24  masses),  4  (16  mosses),  and  7  (35  masses).  The  hull 
motion  is  the  same  in  Cases  3  and  7  but  differs  some  in  Case  4.  The  propeller  motion  (which 
is  really  of  secondary  interest)  is  identical  in  Cases  3  and  7  except  for  a  slight  hump  at 
about  20.  5  cps  in  Case  3.  There  ore  slightly  more  differences  between  Case  4  and  Cases  3 
and  7  in  the  form  of  slightly  sharpened  resonances,  but  these  are  not  large.  It  is  probably 
better  to  use  more  than  16  stations  in  defining  the  hull. 

When  Cases  5,  6,  7,  and  8  are  compared,  it  becomes  clear  that  the  damping  of 
the  propeller  and  the  sprung  masses  is  not  sufficient  to  keep  the  amplitude  of  hull  vibration 
to  reasonable  limits.  However,  the  addition  of  4%  of  hysteresis  damping  is  very  effective. 
10%  and  25%  damping  make  the  submarine  dead. 

Cases  10,  7  end  9  illustrate  the  effects  of  the  thrust  bearing  stiffness  upon  the  hull 
response.  In  Cases  7  and  9  the  effect  of  the  »hrus!  bearing  stiffness  in  the  amplitude  of 
motion  of  the  propeller  is  large  but  its  influence  on  hull  motion  is  small.  However,  Case  10 
shows  that  when  resonant  frequency  of  the  propeller  on  its  thrust  bearing  matches  a  natural 
frequency  of  the  hull,  the  amplitude  of  the  hull  motions  can  be  increased.  When  this  is 

9 

recognized,  it  becomes  clear  that  the  amplitude  of  the  hull  motion  caused  by  relatively 
high  frequency  (>15  cps)  propeller  excitation  is  strengthened  by  propulsion  systems  that 
have  high  axial  resonant  frequencies  (i.e.,  stiff  thrust  blocks).  It  is  probable  that  the 
very  low  natural  frequency  that  can  be  obtained  by  supporting  the  thrust  shoes  on  an  air 
spring,  could  be  very  effective  in  attenuating  the  amplitude  of  hull  vibration  in  the 
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axial  modal. 

Cat#  11  wai  run  primarily  for  estimating  the  amplitude  of  hull  motion  under  sen* 
vice  conditions.  However,  when  compared  with  Cate  7,  it  illustrates  that  a  harmonic 
force  generated  at  the  propeller  will  generally  but  not  always  give  smaller  ship  motions 
than  one  generated  on  the  surface  of  the  hull. 

2.  Vertical  Bending  Vibration 

The  most  conspicuous  influence  upon  the  bending  vibration  is  the  inclusion  of 
the  propulsion  system  as  a  sub-system  in  the  calculations.  This  influence  ik  first  strongly 
shown  even  when  the  excitation  is  on  the  ship  hull  by  comparing  Cases  13  and  14  and  when 
the  excitation  is  transferred  to  the  propeller,  Case  15,  results  in  a  hull  response  that  is 
entirely  different  from  that  obtained  when  the  propulsion  sub-system  is  not  considered. 

A  comparison  among  Cases  15,  16  and  17  in  which  the  hull  is  represented  by 
25,  16  and  35  concentrated  masses  shows  that  the  predominant  propeller  frequency  is  much 
lower  in  Case  15  than  in  Case  16  or  17.  This  is  because  of  the  difficulty  of  tying  the 
propulsion  sub-system  into  the  hull  properly  when  required  to  work  with  regular  station 
spacings. 

In  comparing  Case  16  with  Case  17,  it  will  be  noted  that  although  the  resonances 
occur  at  about  the  same  frequencies,  the  response  is  different  even  at  tow  frequencies. 

This  result  was  not  expected. 

In  comparing  Cases  17  and  18  in  which  the  treatment  of  the  water  inertia  effects 
is  changed  (See  Appendix  C  for  an  explanation  of  this)  it  is  clear  that  no  great  change 
occurs  because  of  the  reduced  inertia  associated  with  the  sinusoidal  approximation  but 
that  the  vibration  does  carry  farther  into  the  hull  at  the  higher  frequencies.  Thus  the 
amplitude  of  vibration  at  the  bow  at  16  cps  is  greater  for  Case  18  than  for  Case  17. 

The  effects  of  hull  damping  are  shown  in  Casas  19  (no  damping)  20  (4%  damping) 
18  (10%  damping)  and  21  (25%  damping).  Referring  to  Case  19  it  is  clear  that  the  hull 
resonances  are  quite  sharp  and  that  the  sprung  masses  and  the  propeller  damping  are 
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Frequency,  cps 

Figure  27  -  Diiplacement  in  Vertical  Bending  for  1  lb.  Excitation  at  Stem,  Ca*e  12A 


Frequency,  cps 

Figure  28-Vertical  Bending  Phase  Angle  Relative  to  Stem  Excitation,  Case  12A 


Frequency,  cps 

Figure  30 -Vertical  Bending  Phase  Angie  Relative  to  Stem  Excitotion,  Case  13 


Figure  31  Displacement  in  Vertical  Banding  for  1  lb.  Excitation  at  Stem,  Cara  14 


32-Vertieol  B«ndi 


Figure  33-  Dhplocement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller,  Com  15 


Frequency,  cps 

Figure  34 Vertical  Bending  Phase  Angie  Relative  to  Propeller  Excitation,  Case  15 


Figure  35-Di*>lacement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller,  Owe  16 
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Figure  37  -Displacement  in  Vertical  Banding  for  1  lb.  Excitation  at  Propel  far.  Cow  17 


Figure  38-  Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  17 


Figure  39- Diiplocement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller,  Cate  18  (Standard) 


^Figure  40- Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  18  (Standard) 


Figure  41  *4)iaplacement  in  Vertical  landing  for  1  lb.  Faroe  of  Excitation  at  Propeller,  Care  19 


iPWHI*  f 


.  Frequency,  cps 

Figure  ^Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  19 


Figure  43-  Displacement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller,  Case  20  (Reduced  Damping) 


Frequency,  cps 

Figure  44-Vertical  Bending  Pha$e  Angle  Relative  to  Propeller  Excitation,  Case  20  (Reduced  Damping) 


Figure  45-  Di placement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller,  Case  21  (Increased  Damping) 
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Frequency,  cps 

Figure  46  Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  21  (Increased  Damping) 


Frequency,  cps 

Figure  47-Displacement  in  Vertical  Bending  for  1  lb.  excitation  at  Propeller,  Case  22 


Frequency,  cps 

Figure  48  Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  22 


Figure  49 -Displacement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller, 


Figure  50- Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  23 


Frequency,  cps 

Figure  51  -Displacement  in  Inertial  Bending  for  1  lb.  Excitation  at  Propeller,  Case  24 


Frequency,  cp* 

Figure 52  Vertical  Bending  Phase  Angle  Relative  to  Propeller  Excitation,  Case  24 
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Figure  53 -Displacement  in  Vortical  Bonding  for  1  lb.  excitation  on  hull.  Com  25 


Frequency,  cp* 

Figure  54-  Vertical  Bending  Phase  Angle  Relative  to  Excitation  on  Hull, 


Figure  56-Verticol  Rending  Phot*  Angle  Relative  to  Excitation  on  Hull 


.Figure  59-  Displacement  in  Vertical  Bending  for  1  lb.  Excitation  at  Propeller, 
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Figure  60 -Vertical  Banding  Phare  Angle  Relative  to  Propeller  Excitation,  Caw  30 
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sufficient  to  hold  them  within  reasonable  bounds  even  without  any  hull  damping.  Surely 
the  effect  of  hull  damping  is  two  fold.  (1)  It  reduces  the  amplitude  of  motion  at  the  point 
o  f  excitation  and  (2)  it  attenuates  the  relative  motion  of  locations  remote  from  the 
source  of  excitation.  The  value  of  damping  which  most  closely  matches  that  of  the 
actual  submarine  will  have  to  be  obtained  by  comparison  with  the  trials  on  the  submarine. 
These  calculations  suggest  that  it  requires  leu  hull  damping  to  match  trial  experience 
than  was  previously  thought.  This  is  a  desirable  result  because  although  experience  in 
other  structures  supports  an  auumption  of  4%  hysteresis- damping,  it  is  very  difficult  to 
explain  the  source  of  10%  or  more  of  hysteresis  damping. 

When  comparing  Cases  22,  18  and  23,  representing  a  10  to  1  change  in  the  stiff¬ 
ness  of  the  stem  tube  bearing  the  change  in  the  hull  response  is  striking.  This  increased 
stiffneu  bearing  (3  times  over  the  standard)  moves  the  propeller  resonance  from  about  11.3  cps 
to  about  1 9. 5  cps  as  might  be  expected.  What  might  not  be  expected  is  the  large  effect 
of  this  shift  upon  the  hull  response.  The  amplitudes  of  motion  at  the  lower  resonances  are 
very  much  reduced  and  those  at  the  higher  resonances  increased  by  this  shift.  The  effects 
of  reducing  the  stiffness  of  the  stem  tube  bearing  are  entirely  consistent  with  the  above 
results.  Although  the  amplitude  of  the  propeller  is  increased  this  occurs  at  lower  speeds 
and  there  is  a  large  decrease  in  the  hull  response  in  the  normal  operating  range.  These 
results  surely  emphasize  the  importance  of  treating  the  propeller  and  shafting  sub-system 
as  precisely  as  possible. 

Case  24  indicates  that  the  hull  is  leu  responsive  to  a  harmonic  moment  (a  moment 
arm  of  more  than  1000  inches  would  be  required  to  give  equivalence)  than  to  a  harmonic 
force.  Cases  25  and  26  as  compared  with  Case  18  indicate  that  except  at  the  propeller 
resonance  a  harmonic  force  at  the  stern  will  generate  a  slightly  higher  hull  response  than 
the  same  force  at  the  propeller.  As  the  location  of  the  harmonic  force  is  moved  forward 
Case  26  as  compared  to  Case  25)  the  response  on  the  hull  falls  off  noticeably. 

Case  27  as  compared  to  Case  18  shows  that  for  the  lower  hull  frequencies  the 
damping  auociated  with  the  sprung  masses  has  a  significant  effect.  It  also  appears  that 
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th.  inelmlon  of  the  .prung  mo»e*  wi»b  tM*  (prlnglng  con  change  the  hull  frequencies  by 
meawreable  amount..  Hoover,  iomo  of  tho  difference.  between  Co...  27  and  18  or. 
difficult  to  under*tand. 

Finally,  Caw  30  a.  compar'd  to  Cce  18,  Indicate,  that  the  effect,  of  including 
the  rotary  inertia  in  the  hull  respome  calculation,  of  a  .ubmarine  are  very  .mall. 
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APPENDIX  A 

PROCEDURE  USED  IN  DEFINING  THE  SHIP  WEIGHT  DISTRIBUTION 

The  basic  information  used  in  determing  the  weight  distribution  of  the  ship  is  given 
in  the  "Detailed  Weight  Summary,  "  SSB(N)  599  dated  7/1/60,  prepared  by  the  Electric 
Boat  Division  of  General  Dynamics  Corporation.  ^  This  weight  summary  is  in  the  form  of  an 
IBM  print-out.  The  summary  breaks  the  ship  weight  into  about  250  classifications.  Under 
each  classification  it  lists  the  plans  covering  the  items  contributing  to  the  weight.  For  each 
plan  it  gives  the  total  weight  represented  by  the  plan  and  the  vertical  and  longitudinal  center 
of  gravity. 

For  many  of  the  weights,  particularly  heavy  concentrated  items,  this  information 
can  be  used  directly  for  estimating  the  longitudinal  distribution.  However,  for  some  impor¬ 
tant  weight  elements,  such  as  the  hull  plating,  the  plans  cover  such  a  large  portion  of  the 
ship's  length  that  it  is  necessary  to  make  independent  calculations  to  define  the  weight  dis¬ 
tribution.  On  many  weight  items,  such  as  electrical  wiring,  piping  and  joiner  work,  the 
distribution  of  the  component  weights  had  to  be  estimated  from  the  compartments  served 
and  the  longitudinal  position  of  the  center  of  gravity. 

From  the  Detailed  Weight  Summary  a  chart  wqs  prepared  in  which  the  weight  in 
each  of  the  several  classifications  was  distributed  throughout  the  ship.  A  typical  section 
of  this  chart  covering  the  portion  of  the  ship  length  between  30  feet  and  80  feet  from  the 
bow  is  shown  in  Figure  A-l.  It  consists  of  distributed  weights  and  concentrated  weights. 

The  distributed  weights  are  added  up  at  the  several  sections  and  a  curve  of  distributed  weight 
prepared.  The  concentrated  weights  are  listed  by  location  and  are  then  distributed  over  a 
two  foot  length.  Table  A-1  presents  a  listing  of  the  concentrated  weights  over  the  interval 
covered  in  Figure  A-1.  Figure  A-2  shows  the  distributed  mass,  the  spread-out  concentrated 
masses,  the  captured  water  and  the  water  ballast  for  the  typical  section  of  the  ship  between 
30  feet  and  80  feet  from  the  bow.  Finally,  all  of  these  curves  are  added  and  smoothed 
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TABLE  A-1  -  SSB(N)  598  -  WEIGHTS  LOCATED  BETWEEN  301  and  80'  FROM 
THE  BOW  THAT  ARE  TREATED  AS  CONCENTRATED  WEIGHTS 


Distance 
from  Bow 
*  ft. 

Weight 

Account 

Number 

Weight 

lbs 

Remarks 

30 

682-»696 

1708 

30.5 

541 -*542 

9015 

Fws.  Escape  Trunk 

33 

344 

5139 

37 

344 

2158 

36 

611-*619 

1348 

36.5 

421 

3358 

41.6 

417 

791 

42 

344 

1152  1 

42 

353 

318 

5241 

42 

421 

1297 

42 

61 1-*619 

2474 

42.5 

433 

223 

44.8 

417 

244 

50 

421 

25 

50.5 

611-*619 

1348 

53.4 

537 

159 

364 

53.4 

682-*696 

205 

54 

543 -*545 

155  ^ 

■ 

54 

641 -*644 

520 

l  12,175 

54 

648 -*649 

11,500  J 

Oxy.  Tanks  and  Piping 

56 

344 

83 

56 

422 

151 

56 

627-*  634 

522 

3693 

56 

451 -*462 

2937, 

56.5 

541 -*542 

768 

56.5 

543 ‘♦544 

155 
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TABLE  A- 1  -  (continued) 


Distance 
from  Bow 
ft. 

Weight 

Account 

Number 

Weight 

lbs 

Remarks 

im 

415 

671  ) 

733 -#738 

176  }  1007 

57.5 

761 -*767 

160  J 

58.6 

331 

54,332 

Trans.  Bkhd. 

59 

991-*999 

2457 

61 

641 -*644 

1984 

74 

561 

12,940 

Lead  Ballast 

74.5 

353 

1383 

75 

421 

936 

A-ll 


Figure  A- 2  -  Section  of  Weight  Distribution  Curve 
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slightly  and  the  resulting  curve  used  to  represent  the  weight  distribution  on  the  ship.  This 
curve  is  given  in  Figure  .3/  .  For  a  check  this  curve  is  integrated  and,  with  captive  water 
deducted,  compared  with  the  submerged  displacement  of  the  ship. 

With  the  general  process  explained,  we  will  now  discuss  some  of  the  details.of 
the  process  of  distributing  the  weights  along  the  hull.  For  some  of  the  main  categories  of 
hull  weight  it  is  necessary  to  compute  the  weight  distribution  along  the  hull.  This  is  not 
particularly  laborious  in  the  case  of  the  inner  and  outer  hull  shell  plating,  since  the  areas 
must  be  computed  in  the  stiffness  determinat:on.  Lead  ballast  is  a  very  significant  item  of 
weight  and  can  be  concentrated  within  a  small  length.  It  is  therefore  well  to  define  this 
weight  closely  from  the  plans.  Also  the  water  ballast  is  a  very  significant  weight  item. 

Normally,  this  is  easily  defined  from  the  weight  in  a  particular  tank  but  in  the  case  of  the 
GEORGE  WASHINGTON  it  was  necessary  to  study  the  distribution  of  weight  in  the  forward 
ballast  tanks.  Items  like  watertight  flats  are  usually  broken  down  into  separate  compartments 
and  can  be  distributed  quite  easily  by  brief  reference  to  the  plans  to  ascertain  the  location, 
extent  and  shape.  The  framing  is  broken  down  in  the  Detailed  Weight  Summary  into  groups 
of  four  or  five  frames  and  the  weights  of  these  are  quite  easily  distributed  over  their  longi¬ 
tudinal  span. 

The  distribution  of  the  weights  under  some  classifications  such  as  furniture,  wiring,  and 
piping  are  more  difficult  to  make.  Where  it  is  clear  from  the  plan  title  that  the  weights 
represented  by  the  plan  fall  within  a  certain  span  of  length,  they  can  be  distributed  either 
uniformly  or  by  some  simple  distribution  over  the  length  of  the  span.  For  convenience  in 
determining  weight  distribution  for  such  weight  categories  as  interior  painting,  electrical 
wiring,  telephone  wiring,  etc.,  the  interior  of  the  hull  was  considered  to  be  divided  in 
ten  sections  as  indicated  in  Table  A- 2.  Using  these  sections  it  is  possible  to  estimate  a  weight 
factor  that  will  indicate  what  proportion  of  the  weight  will  fall  within  each  section.  The 
application  of  this  process  to  the  interior  paint,  account  261,  which  weighs  28,  207  pounds 
is  illustrated  in  Table  A-3.  The  application  of  the  process  to  the  distribution  of  the  weight 
of  power  wiring  through  the  ship  (weight  about  45,000  lbs)  is  illustrated  in  Table  A-4.  The 


Hull  Weight  SSp(N)  598  GEORGE  WASHINGTON  t 


CONESCO,  INC.,  consultants  In  engineering  science 


TABLE  A- 2  -  HULL  DIVISIONS  FOR  ESTIMATING  WEIGHT  DISTRIBUTION 


Section 

— 

Compartment 

Fr.  Nos. 

Distances  from  F.  P. 

■■9 

0 

Bow 

0  — ♦  13 

0  -♦  25. 25 

25.25 

1 

Torpedo  Rm.  Fwd. 

-»  20 

-»  40.83 

15.58 

2 

Torpedo  Rm.  Aft. 

—  26 

-*  57.50 

16.67 

3 

Living  Qtn. 

—  M3 

->  92.75 

35.25 

4 

C.  O.C. 

-»  Ml  7 

-*  1 32. 92 

40.17 

5 

Missile  Rm. 

— ♦  M44 

-» 208. 92 

76.0 

6 

Missile  Compt.  Aft. 

-*  44 

-*  233. 67 

24.75 

7 

Reactor  Compt. 

-»  52 

-*  253. 67 

20.00 

8 

Aux.  Mach.  Space 

-»  64 

-t  285. 67 

32.0 

Engine  Rm.  Fwd. 

_  77 

->320. 17 

34.50 

H 

Engine  Rm.  Aft. 

—  92 

~>  353. 67 

33.50 

Stern 

—  102 

-►373. 50 

19.83 
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TABLE  A- 3  -  DISTRIBUTION  OF  WEIGHT  OF  INTBUOR  POINTS 


1— 

longth 

Product 

wt/ft 

1 

15.6 

0.7 

10.9 

2 

16.7 

1 

16.7 

3 

35 

3 

105 

4 

40 

3 

120 

145 

5 

76 

2 

152 

97 

6 

25 

1 

25 

48 

7 

20 

1 

20 

48 

8 

32 

1.5 

48 

73 

9 

34 

1.5 

51 

73 

10 

34 

1 

34 

48 

583 

wt/ft  =.  x  -  48.4  weighting  factor 
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TABLE  A-4  -  DISTRIBUTION  OF  WEIGHT  OF  POWER  WIRING 


Assume: 

a  1  cable  from  compartment  2  to  compartment  9 

b  1  cable  from  compartment  3  to  compartment  9 

c  2  cables  from  compartment  4  to  compartment  9 

d  1  cable  from  compartment  5  to  compartment  9 

9  2  cables  from  compartment  6  to  compartment  9 

f  1  cable  from  compartment  7  to  comportment  9 

g  2  cables  from  compartment  8  to  compartment  9 

h  3  cables  in  compartment  9 

i  3  cables  in  compartment  10  to  compartment  9 

wt/ft 

compartment 

ab  c  d  efgh  i 

Total 

25 

2 

17 

17 

50 

3 

35  35 

70 

99 

4 

40  40  80 

160 

124 

5 

76  76  1  52  76 

380 

148 

6 

25  25  50  25  50 

175 

198 

7 

20  20  40  20  40  20 

160 

248 

8 

32  32  64  32  64  32  64 

320 

320 

9 

34  34  68  34  68  34  68  102  100 

440 

73 

10 

100 

100 

1822 

.  .  I  oral  tv,  wv 

wt/ft  in  compartment  =  6rqn<fTo>-qr  x  gf  cqffyqrf— y 
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application  to  othor  categories  such  at  tho  announcing ‘system,  tho  telephone  system,  otc.,  - 
can  bo  understood. 

Equipment  weights  which  exceed  2000  pounds  and  are  shock  mounted  are  treated 
as  sprung  masses  (See  Appendix  G).  The  fairwater  planes,  the  stabilizers  and  the  stem 
planes  are  especially  studied  as  elastic  structures  in  vertical  vibration  but  are  considered 
as  rigid  masses  in  the  axial  direction.  The  rudders  are  considered  to  be  rigidly  connected 
in  both  vertical  and  longitudinal  vibration.  Properly,  the  stabilizers  and  control  surfaces 

|  should  have  been  incorporated  in  the  calculation  program  as  elastic  sub-systems  but  since 

I 

the  program  was  not  developed  when  the  hull  was  defined,  this  was  not  done. 

I  The  "captured"  water,  water  in  the  free  flooding  sections  of  the  ship,  is  treated 

| 

as  part  of  the  mass  of  the  ship.  Consider,  for  example,  the  water  between  the  superstructure 
and  the  hull.  Beneath  the  superstructure  the  water  has  a  relatively  free  flawing  space. 
However,  the  air  vent  in  the  superstructure  deck  is  only  a  very  small  percentage  (leu  than 
5%)  of  the  superstructure  deck  area.  Since  it  normally  passes  air  this  area  is  adequate  for 
venting  but  water  flawing  through  it  would  have  to  pass  at  more  than  20  times  the  mean 
velocity  in  the  tank.  As  a  consequence  for  any  vibration  when  the  submarine  is  under  water 
the  water  cannot  flow  through  the  superstructure  vents  in  sufficient  quantity  to  allow  any 
significant  free  flew  and  the  water  is  indeed  "captured"  and  forced  to  move  with  the  ship. 
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APPENDIX  B 


VIRTUAL  INERTIA  OF  ENTRAINED  WATER  IN  LONGITUDINAL  MODES 


1.  Determination  of  an  Ellipsoid  to  Represert  the  Hull 

The  distribution  of  the  inertia  of  the  entrained  water  on  the  SSB(N)  598  is 
assumed  to  be  the  same  as  that  of  the  ellipsoid  of  revolution  which  most  closely  fits  the 
hull.  The  first  step  in  determining  the  inertia  of  the  entrained  water  therefore  is  to  find 
the  dimensions  of  the  ellipsoid  which  most  closely  marches  the  hull.  The  process  that  is 
used  is  the  "Method  of  Least  Squares"  and  is  described  in  Reference  8. 

Although  the  axes  of  the  portions  of  the  submarine  that  are  surfaces  of  revolution 
are  offset  in  the  vertical  direction  and  some  of  the  coss-sections  are  not  circular,  it  was 
assumed  that  the  prolate  ellipsoid  which  matches  the  breadth  of  the  submarine  most  closely 
would  adequately  represent  the  submarine.  The  origin  of  the  axes  of  coordinates  was  taken 
at  the  centerline  at  the  bow  and  the  x  direction  was  taken  along  the  centerline  of  the  bow^ 
and  stem  sections.  The  equation  of  the  ellipse  that  rotates  into  the  ellipsoid  is  ^ 
where  a,  b  and  c  are  to  be  determined  by  the  method  of  least  squares  to  give  the  a 
ellipsoid  that  most  closely  approximates  the  submarine. 

Assume  for  a  first  approximation  that  a  -  b  =■  190'  and  c  =  16. 5'. 

From  the  equation  of  the  ellipse 


y  =  c 


’/ 2 


Expand  this  value  of  y  in  a  Taylor's  series. 


y  -y, 


a,b,c 


a,  b,  c 


O,  b,  c  a,  b,  c 


Ac  +: 
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where: 


Letting  s  equal  the  sum  of  the  squares  of  the  differences,  Vj ,  between  the  actual  value 
of  y,  y. ,  and  the  value  computed  by  the  equation  for  y. 


The  values  of  Aa,  Ab  and  Ac  which  will  minimize  s  are  found  by  setting 


ds  8s  8s  n 

l^A a)  *  "8(Ab)  “  *9(Ac)  =  t 
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This  /(aids  the  normal  aquations: 


n  l"  .  .  21  ‘  ]//2  .  ,.2 

■Vi"*?  *Tl- 

L  a  a 


n 

21  v.  c 


fca  =  o 


■  [,-^2J  1/2  =0 

After  substituting  the  value  of  v. ,  these  equations  reduce  to: 

vi  ut2  A  c  a  ?  f,  fr-b)*l  \ 

c  I  (x-b)  +  -^  Aa  I  1  -  5 — jS-  (x-b) 

a  L  a  _ 


3  n  o 

(x-bj  +  Ac  Z  (x-b) 


(x-b)  =  0 


?4-¥\ 


cl  (x-b)  +  -3 
a 


♦  AbZ  1  -  (x-b)2  +  Ac  I  (x-b) 
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To  find  the  values  of  Aa,  Ab  and  Ac,  tablet  are  formed  for  the  quantities 

n  n  o 

I  (*-b),  I  (x-b) 


e<2o) 


B(3a) 
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Values  of  y.  are  picked  from  the  graph  of  the  half  breadths  at  the  15*  -  9-7/8"  water  line 
as  plotted  in  Figure  B-1  for  six  values  of  x. 

The  results  are  given  in  Table  B-1. 

Substituting  the  values  from  Table  B-1  into  the  equations  B(1a),  B(2a)  and  B(3a)  yields  the 
following: 


17,380  Aa  -  5690  Ab  +  88,700  Ac  =146,000  B(lb) 

0.2138  Aa  -  0.00456  Ab  +  3.543  Ac  =  2.0  B(2b) 

29.95  Aa  -  132  Ab  +  10  Ac  =  2152  B(3b) 

From  a  simultaneous  solution  of  these: 

Aa  =  0.517'  a  improved  =  190‘  +  0.517"  =  190.52' 

Ab  =-16.12'  b  improved  =  190'  -  16.12'  =  173.88’ 

Ac  =  0.515*  c  improved  =16.50'  +  0.515'  =  17.015' 


Therefore  the  ellipsoid  that  most  closely  matches  the  submarine  is; 

(x-173.  88) 2  +  (y2  +  z2)  =  } 

(190.  52) 2  (17.015)2 

The  curve  of  this  quantity  is  plotted  in  Figures  B-1  and  B-2.  It  is  obvious  that  although 
this  ellipsoid  matches  the  ships  hull  as  closely  as  is  possible  by  any  single  ellipsoid,  there 
is  still  a  considerable  discrepancy  at  the  ends. 

Because  of  the  complexity  of  dealing  with  a  closer  simulation  of  the  ship,  the 
simple  ellipsoid  has  been  considered  to  be  an  adequate  representation  at  the  present  time. 
The  methods  of  Reference  9  would  have  to  be  developed  further  before  they  would  be 
applicable  to  the  problem  of  a  shape  differing  from  an  ellipsoid. 


TABLE  8-1  -  MEAN  SQUARE  DETERMINATION  OF  ELLIPSOID 
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E 

<0  *0  'O  'O  *o  *o 

o  o  o  o  o  o 

5  5  5 

8  fc  T  8  8 

-  »0  N  ^  (M  CS 

V  csT 

1 

B 

°o  *0  '°o  *°o 

"x  X  ~X  X  7  7 

«  2  S  ?  8  2 

•*  .  <5  cm  .  . 

CM  V  CM  sr 

1  1  o  +  — 

• 

s 

? 

cs 

i 

E 

144,500 

31,780 

1,670 

3,990 

18,200 

88,400 

s 

rc  ■ 

“?  <N 
^  C 

0.20 

0.  531 

0.956 

0. 902 

0.664 

0.290 

s 

co 

CN 

J? 

1 

A 

28.900 

16.900 

1,600 

3,600 

12, 100 

25,600 

8' 

S 

Jc, 

R  8  ?  S  2  8 
7  7  ’*77 

_ 

1? 

O  O  »n  O'  i— 

o’  in  o’  n  ri  is 

X 

8  3  2  §  §  § 

CM 

JO 

A 

CM 

N 

T 

1  T 

1 _ l 

X 

646,000 

347,000 

27,000 

58,300 

191,300 

338,000 

8 

% 

2T 

i 

3S. 

<* 

rJT. - - 

1 _ 1 

X 

-3800 

-2670 

-675 

+973 

1,740 

2,115 

N 

W 

cm' 

1 

> 

1 - ! 

1 _ 1 

X 

4.47 

10.93 

16.15 

14.62 

10. 51 

3. 82 

8 

• 

8 

1  04  i 

\ 

fh> 

i 

1 _ !I_l 

.447 

.730 

.978 

.950 

.815 

.538 

X 

X 

10.0 

15.0 

16.5 

15.4 

12.9 

7.1 

8  9  2  §  8  8 
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2.  Analytical  Representation  of  the  Vibration  Pattern 

It  is  also  necessary  in  the  considerations  of  equivalent  mass  to  work  with  a 
pattern  of  vibration.  For  this  purpose  the  nodal  pattern  at  the  natural  frequency  of  10.40  cps, 
as  computed  for  the  mass  and  stiffness  reported  in  Reference  3  was  used.  This  nodal  pattern 
is  plotted  in  Figure  B-3.  For  convenience  in  the  computations  a  deflection  curve  of  the 
form  6  =a}  (173.  88  -  x)  +  a3  (173.  88  -  x)3 

*  OjX'  +  a3x'3 

was  assumed  and  the  values  of  a^  and  a3  to  match  the  experimental  curve  determined  by  the 
method  of  Least  Squares.  The  resulting  analytical  expression  for  the  deflection  is: 

6=4.11  x  10"6-  11.08  x  10”1 2  (x1)3 

This  is  plotted  on  Figure  B-3  to  indicate  the  adequacy  of  the  match. 

3.  The  Distribution  of  the  Entrained  Water  Inertia  on  the  Ellipsoid 

We  are  now  concerned  with  the  water  ineirtia  of  an  ellipsoid  of 
major  axis,  a^  -  190.52' 
minor  axes,  b^  and  c^  -  17.015' 

•6  “12  3 

vibrating  in  a  longitudinal  mode  such  that  6  —  4. 1 1  x  10  x  -  11.08  x  10  x 
where  x  is  measured  from  the  center  of  the  ellipsoid. 

The  velocity  in  the  longitudinal  direction  is  -  jwS 
Express  this  as : 

u  =  Ax  -  Bx3 

e 

Because  of  the  poisson  ratio  effects  there  will  be  a  rodial  velocity  to  the  hull  associated 
with  the  longitudinal  velocity. 

36 

The  longitudinal  strain  =  — 


B-9 


figure  B-3  -  Representation  of  Normal  Mode  Deflection  by  a  Functional  Relationship 
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fc 


d& 

Corresponding  to  this  longitudinal  strain  there  will  be  a  poisson  strain  ?  ££ » 

where  o  is  a  ratio  that  depends  upon  the  poisson  ratio  of  the  material  of  the  submarine  and 

—  86 

the  ratio  of  frame  area  to  shell  area.  Hence  the  poisson  velocity  *  y  (or  z)  a  ^  and 

v  =  ay  (A  -  3  Bx^) 

© 

w  =  erz  (A  -  3  Bx^) 

© 

The  methods  used  for  determining  the  water  inertia  were  originally  developed  by 
Lamb^®  and  are  well  presented  in  Reference  1 1.  The  methods  involve  the  expression  of  the 
properties  of  the  ellipsoid  and  the  vibration  in  terms  of  ellipsoidal  coordinates.  In  this 
manner  a  potential  function  can  be  obtained.  In  most  works  with  virtual  inertia  it  is 
possible  through  energy  considerations  to  obtain  an  equivalent  mass  for  the  water  vibrating 
with  the  ellipsoid.  However,  in  this  case  the  effects  of  the  entrained  water  are  distributed 
over  the  full  length  of  the  ellipsoid  and  in  order  to  find  the  effects  of  this  entrained  water 
at  any  section  of  the  hull  it  is  necessary  to  compute  at  each  longitudinal  location  the  com¬ 
ponent  of  the  force  associated  with  the  dynamic  pressures  between  the  hull  and  the  woter 
that  result  from  the  modal  pattern.  This  force,  when  related  to  the  amplitude  of  vibration 
at  this  location,  can  be  used  to  determine  the  equivalent  moss  of  the  water  at  this  location. 

The  discussion  of  the  geometry  and  the  kinematics  of  a  vibrating  ellipsoid  given 
in  Reference  11  is  directly  applicable  to  this  problem  and  therefore  will  not  be  repeated. 
Because  the  concern  of  this  study  is  axial  vibration  rather  than  bending  vibrations  the 
applications  of  the  general  theory  differ. 

Express  the  velocities  in  ellipsoidal  coordinates 

*  a  3 
u  =  Ax  -  Bx 

e 

=  Aku?  -B(ku?  )3 
o  u 
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v#  =  cry  (A  -  3  Bx2) 

=  ak  (1  -  u2)  1/2  (?  /  '  1)  1/2  co*  6  (a  -  3B  (kuC  o)2] 
w#  =  ok  (1  -  u2)  ^2  (?o2  -  1)  ^  sin  ej^A  -  3B  (ku?Q)2] 

Since  these  velocities  are  those  in  the  surface  of  the  ellipsoid  this  indicated  by  the 
ellipsoid  representing  the  hull. 

2  'A  2  2 

FU9)  =  yu#  +  (1-|T)  ?Q  (?Q-1)  (ve  cos  6+we  sin0) 

-i*  [*•>»*«.- 

9  Vo  5  Vi  r  2 

+  0  -  H  )  2  ?0ok(1-/)  I^A  -  3B  (k^? q)  1 

=  p  AkV?o-B(kH?o)3  +  ok?o(l  -y2)  A-3B(kH?o)2 

=  A  kf2?o  +®k?o  (1  -p2)  -  B  (ku?Q)3  +3ok?o(l -p2)  kp?o2 


=  A  ku2?o  (1  -  0) +ok?Q  -  B  (kp?o)3  ^i(l  -  3o)  +  ^  o  (k^?Q)3 


Applying  the  condition  that  k  F  (^,  6)  =  ,  to  the  general  solution  for  velocity  potential 

in  ellipsoidal  coordinates,  i.e.. 


*  =  I  I  (B  *  sin  s9  +  C*  cos  s0)  P*  (u)  Q*  (5  ) 

n  _q  j  ^  n  n  n  *  n 
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Uadi  to 


n 

I  £ 

n  =0  i  <*0 


(B*  »ln  >6  +  C*  coi  i9)  Pn*  (H)  Q*  <g  - kF  0) 


Sine*  the  F  (^,  0)  It  Independent  of  0  thii  condition  ii  only  met  where  i*( 

Therefore  B*  falls  out  of  the  solution 
n 


C  becomes  C 
n  n 


P  becomes  P 
n  n 

~S 


Hence 


Q  becomes  Q  and  the  equation  becomes, 


ZCnPn  ¥  Qn«o>=kFV'e> 
n  =U 


c =  kx2w  f  F  (u,  0)  P  ( J  d^i 

4*Qn(?0)  A  r 

no  “l 


=  2n+1 
2Q  «  ) 

n  o 


i 

Ak?o(l-a)  J  p2Pn(^)du 
-1 


I  I 

A’k!o  f  W  Jr  Bk  V  ,,-3<’)  f  f<Pr 

-1  'I 


1 

f 


■1 
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And  using  Table  1  of  Reference  1 1 


C  = 
o 


2QA> 


[Ak? 


(2  +  4o)  -  Bk3 5 Q3  (3  +3C) 


C1  =  ° 


C„  = 


2GW 


k  AkS0(w)"Bk3?<>3('3r  +  '75'  a) 


c_  =  0 


C4  " 


2Q4  (?  0) 


D|3  u  3  i.  n  s  I 
-Bk  (1-3.)  3,3- 


C  =  0,  n  i  5 
n 


Therefore: 


<P  = 


Q  (?J 

O  O 


k  (1  +  2a) 


Ak?  Bk3?  3 
o  o 


po  (fj  Q0  (? ) 


<vy 


k  §  Ak?o  (l-a)  -^Bk3?o3  (2+a) 


P2^i)Q2(?) 


srr.  k*V<i-a>)  p4v  q4«> 


«y 
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When  the  potential  function  of  the  fluid  It  determined,  the  usual  practice  for  obtaining 
the  mau  of  the  entrained  water  it  to  compute  the  kinetic  energy  of  the  fluid  and  by  re-* 
lating  this  to  the  amplitude  of  motion  to  determine  the  equivalent  matt.  However,  for 
this  problem  the  distribution  of  the  equivalent  mats  along  the  length  of  the  hull  it  de¬ 
sired,  and  the  amplitude  of  vibration  is  not  a  constant  value  but  a  function  of  the  longi¬ 
tudinal  position  along  the  hull. 

To  avoid  these  difficulties  the  distribution  of  the  equivalent  mass  is  obtained  by 
computing  the  pressure  on  the  ellipsoid  generated  by  the  vibratory  motion,  multi  plying 
this  pressure  by  the  projected  area  normal  to  the  longitudinal  axis  for  each  unit  of  axial 
length.  The  equivalent  mass  per  unit  length  at  a  given  axial  position  is  then  the  force  per 
unit  of  length  divided  by  the  product  of  the  amplitude  of  vibration  at  the  location  and  the 
square  of  the  circular  frequency. 

The  pressure  associated  with  a  velocity  potential  function  is; 


p  =p 


d<p 

IT 


+ 


\  +  f(t) 


In  this  relation  q  is  the  square  of  the  fluid  velocity,  which  for  the  small  motions  is  negligible. 
f(t)  is  an  external  pressure  that  does  not  exist. 

Thus  for  this  problem 


d<J>  .  . 

p  =  P37  =  Pl^ 


since  the  ellipsoid  is  vibrating  in  simple  harmonic  motion. 
Since 

2  '2 

a=k*0  c  =  <C-n 

0  0 
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JL 

jup 


Q  (?  ) 

o  o' 


<¥'o> 


Va2-  c2  (1  +  2c)  ^ 


Bo' 

T 


Po^ao«o) 


7a2 -c2  |  Aa  (1-a)  -  ^  Ba3j<2+c)]  P2(H)  Qj  (?( 


8  - — 1 V  a2-  c2  Ba3  (1- 

Q4<V 


The  axiol  force  =  p  x  tr 


0-3o)  P4(H)Q4(?0) 

2  ,  .  dy  ,.2] 


y  -  (y  +  g  xi)' 


dy 


=  -px»x2y^  xl 


Since  for  an  ellipsoid 


dy  xc 

^d x  "  "  ~T 
a 


The  axial  force  =  p  (2w  -y  x  )  lb  per  ft. 


2wp  -yX 

The  entrained  inertia  =  a 


u26 

The  mathematical  expression  for  the  deflection  curve  which  was  derived  earlier  is 

.  3 

6  =  a,  x-a3x 
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or  in  term*  of  velocity  coefficients; 


S 


x 


3 


Tho  entrained  inortia  tharafora 


2ap  c_ 


A-Bx' 


P 

A  _  Ba2 

C2 

(1  +  2a) 

o 

o 

-a, 

T  7* 

L  J 

Po  <H>  Qo  «o’> 


<W 


\  A  (1  -  a)  -  y  Bo2  (  2  +  a  ) 


P2  <H>  Q2('o> 


8  1 
15 


Bo  (1  -  3a)  P4  (v)  Q4(?o) 


where; 


x  o  2  .22 

H  =  a  '  B*  =Ba  »* 


For  the  submarine  SSB(N)  598  a  =  190.52’ 

c  =  17.015’ 


?  2  =  (190. 52) 

0  (190.52)2- (17.01 5)2 


r  =  1.004012 


=  =,'W804008 

r  3  =  1.012084 
»  o 

r  4  =  1.0161448 

»  Q 
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Qo<n 


i ,  ?  +1 

2lofl  f-rj 


Q2(?) 


\  (3?2-  1)  Qo(0  -  3/2  ? 


Q4(?)  =  3  (35?  4  -  30?  2  +  3)  Qo(?)-  ^?3  +  H? 


Q0(?) 


1  i  ~  1  ff  -  1)  1  -  ft  »  1)  * 
1  <?.  ,)2 


?2- 1 


Q2(n 


=  \  (6?  )  Q0ft)  +  5  (3? 2  -  1)  -3/2 


=  3?  Q  (5)-  %  — 
°  r2-  1 


Q/) 


-Jo».«»-«k,q0«-;  (35;4-f2t3> 

r  - 1 


105  „ 2  ,  55 
T  5  +  W 


Q4ft)  = 


35  _ 3  15 _ 

T  »  ' 


Qo<?>  -i 


35? 4  -  30?  2  +  3 
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♦  lost4-  T05C2  -  “t2  ♦  “ 


t  -1 


^t)oo(t)-  '  (3St4-!^  ♦¥> 


t2-. 


-f  Pt3-3t)  00(t)-j  (I05t4-  use2  +  16) 


t2-' 


Q„  «„)  -  5  slw  -  0.6951  -<1.3893  -  6. 9078)] 


=  -j  (6.2136)  =  3.1068 


q2(?o)  *  j  (3  x  1.00804-  1)  x  3.  1068  -  3/2  (1.00401) 


=  1.6382 


q4(Co)  *  3  [35  x  (1.00804)2  -  30  (1.00804)  +  3 


35 


3.1068-^(1.00804) 


=  3.23257 

-  4.42787  =  1.10556 
+  2.30086 


Q  (?)=-  — 

o  o  « 2 


?‘-1 


.0553555 


=  -  124.377 
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Q2  (?q)  =  3x1.004012  x  3.1068-  124.377(3x  1.00804-2) 

=  9.35768-  127.377 
=  -  118.019 

Q4(?o)  =  |  (7x1.012084  -  3  x  1.004012)  3. 1068 

-  j  (105  x  1.0161448-  115  x  1.00804008+  16)  ^864008 

=  31.6315-  140.3509 

=  -  108.7194 

b2  =  (17. 015)2=.  289.  51  /a2  -  c2  =  V36,008.36 

=  189.76 


a  =  approx.  0. 20 
The  entrained  inertia  = 


2»p  x  289. 51  x 

190.52 


(1.40)  x 


3. 1068 

-7547177 


A 

7 


.  1.6382  f  2(0.8) 
-118.019  [  T~ 


(2.2)  Ba 


8 


'  15 


1.10556 

^775775 


(1  -  0. 6)  Ba2  P4  ((i) 


1 


1 

A-Ba 


TT\ 

H 


Po<H> 


B-20 


CONESCO,  INC./  comultanti  In  engineering  science 


1.99 


lb  uc1 
ft4 


11,08x10  ' 
4. 1 1  x  10”6 


x  (190. 52)2 


=  0. 0976  for  2nd  mod*  of  SSB(N)  598 


The  entrained  inertia  then  becomes; 


3600  0. 01093  Pq  (ja)  +  0. 00655  ?2  (^)  -  0. 0000976  ?A  ^  J 


1 


1-0.0976* 


-39.4  [poV 


+  0.  599  P„  (u)  -  0.00892  P.,H/  - , 

r  4  r  1  1-0.0976/ 


’4  V] 


1 


The  values  of  this  entrained  inertia  are  computed  in  Table  B-ll  and  plotted  in  Figure  B-4. 


4.  The  Distribution  of  Entrained  Water  Inertia  on  the  Actual  Hull 

Since  the  shape  of  the  hull  differs  from  the  ellipsoid  it  can  be  expected  that  the 
inertia  of  the  entrained  water  will  differ  from  that  on  the  idealized  ellipsoid.  To  investi¬ 
gate  haw  much  difference  this  might  be,  assume  that  the  pressure  generated  on  the  hull  at 
any  position  along  the  longitudinal  axis  is  the  same  for  the  actual  submarine  as  for  the 
ellipsoid.  When  this  is  done  the  distribution  of  water  inertia  is  quite  different  as  is  shown 
in  Figure  B-5.  It  will  be  noted  that  the  blunt  bow  and  the  discontinuities  at  the  super¬ 
structure  tend  to  give  large  values  of  water  inertia  and  the  pointed  stern  tends  to  decrease 
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TABLE  B-ll  -  CALCULATION  OF  LOCAL  VALUES  OF  ENTRAINED  MASS 


a  r 

p2<m> 

-.  00892P^(fj) 

IB 

— 

■ 

K 00976 

yfli 

0 

-.50 

-.300 

0.375 

-.0033 

0.697 

1 

27.43 

1 

1 

0.599 

1 

-.00892 

1.590 

.9902 

63.30 

0.5 

-.125 

-.0749 

-.28906 

.0026 

0.9277 

.9976 

36.60 

0.75 

. 34375 

.2057 

-.35010 

.0031 

1.2088 

.9945 

47.9 

0.  25 

-.40625 

-. 2438 

.15771 

-.0014 

.7548 

.9994 

29.78 

Total  Entrained  Water  Mass 


H 

Mass/ft 

SM 

Product 

lam 

27.43 

1 

27.43 

0. 25 

29.78 

4 

119.12 

0.5 

36.60 

2 

73.20 

0.75 

47.90 

4 

191.60 

1.0 

63.30 

1 

63.30 

Total 

474.65 

Mass  =  474. 65  x  x  2  =  1 5,  380  lbs.  This  is  identical  with  the  gross  value  of 

entrained  mass  used  in  the  longitudinal  vibration  calculations  in  Reference  3. 
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the  water  Inertia.  Because  of  the  large  difference  between  the  water  inertia  given  by 
Figure  B-5  and  that  of  an  idealized  ellipsoid  shown  in  Figure  B-4,  it  would  be  highly 
desirable  to  develop  an  accurate  method  of  computing  the  water  inertia  on  an  irregular 
body  but  the  procedures  for  doing  this  are  not  developed  and  are  outside  of  the  scope  of 
this  project. 

For  the  longitudinal  vibration  studies  on  the  SSB(N)  598  the  values  of  water  inertia 
given  in  Figure  B-5  are  used.  It  should  be  noted  that  in  the  longitudinal  modes  the  total 
mass  of  the  entrained  water  is  only  15, 000  pounds  (for  either  the  idealized  ellipsoid  or  the 
actual  hull  shape).  This  is  about  1%  of  the  submerged  weight  of  the  hull.  The  small  value 
of  this  entrained  mass  contrasts  strongly  with  *he  large  values  of  entrained  mass  in  the 
bending  modes. 

5  A  Brief  Consideration  of  the  Effects  of  Fluid  Elasticity  upon  the  Entrained  Mass  and 
Energy  Dissipation  of  a  Prolate  Ellipsoid  Vibrating  Longitudinally 

When  a  ship  is  vibrating  in  its  bending  modes  it  has  been  shown  that  there  is  no 
significant  energy  dissipated  in  the  water.  This  can  be  explained  in  terms  of  the  length  of 
a  wave  length  in  the  ship  at  a  given  vibration  frequency  as  compared  with  that  of  sound.in 
the  water  at  '•he  same  frequency.  It  can  also  be  shown  that  the  entrained  water  carried  by 
the  ship  in  bending  vibration  can  be  computed  accurately  without  considering  the  com¬ 
pressibility  effects  of  the  water.  However,  when  a  submarine  is  vibrating  axially  it  is  known 
that  vibratory  energy  is  fed  into  the  water  -  explained  in  terms  of  the  wave  length  of  the 
axial  vibration  being  much  closer  to  the  wave  length  of  the  sound  in  the  water.  For  this 
reason  it  is  desirable  to  have  a  solution  for  the  entrained  mass  and  for  the  energy  dissipation 
of  an  ellipsoid  vibrating  longitudinally. 

The  solution  for  the  vibration  of  a  prolate  ellipsoid  in  a  compressible  fluid  is 
similar  to  that  for  the  incompressible  fluid  but  somewhat  more  complicated.  Where  the 
fluid  is  incompressible  the  equation  for  fluid  potential  is  v  L ^  =0  but  when  the  fluid  is 
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2  | 

compressible  the  scalar  wave  equation  v  41  =  — *-  — v  «|»  applies. 

C  8t; 


For  simple  harmonic  motion 


*  (x,y,  z,  t)  =  4j  (x,  y,  z)  e 


-iuT 


and  the  equation  simplifies  to  the  scalar  Helmholtz  equation 


V  2 «l<  +  p2  <\>  =  0,  £  =  U/C 


BW 


For  the  submarine  problem,  it  is  desirable  to  express  this  equation  in  ellipsoidal 
coordinates  (u,?  ,  9).  These  are  related  to  the  rectangular  coordinates  by: 


x  =  k|J? 

9  l/o  9  Vo 

y  =  k  (1  -f)  1)  *  cos  0 

2  ^2  2  ^2 

z=k(l-p)  (?  -  1)  sinG 


1  £  ?  ^ 


-1  -  ii  ^1 


0  s=  e^2w 


B(5) 


foci  are  at  x  =  +  ^  k 


Then  the  Helmholtz  equation  takes  the  form: 


_8_ 

8V 


2  2.  2,2  . 

-  (?  ”  y  )  f  k  4* 


B(6) 
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This  may  alio  be  writtan  In  the  form: 


Whan  we  let.^i  -  <p  (0)  S  W  J  (?) 

The  equation  separates  to  the  ordinary  differential  equations 

<r+m2<*>-0  B(7a) 


d 

W 


(A-h252 


4 


J  -  0 


B(7b) 


d 

¥ 


L(,V)|-J  +^-hV-  fr )  5=0 

r  i-p 


B(7c) 


where  h  =  and  m  and,  A  are  separation  constants. 


The  solutions  of  B(7a)  are  <p( 9)  —  cos  m@ 

sin  m0 

Equations  B(7b)  and  B(7c)  are  identical  and  so  the  angular  function  $  and  the  radial  function 

J  are  solutions  of  the  same  equation  over  different  ranges  of  the  variable.  These  equations 

can  be  carried  too  solution  through  the  use  of  spheroidal  wave  functions. 

12,13,  14, 

Studies  of  the  radiation  of  noise  from  submarine  hulls  show  that  for  the 
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frequency  range  of  propeller  generated  excitations, the  influence  of  compressibility  of  the 
water  is  small  and  will  tend  to  decrease  the  amount  of  virtual  mass.  Since  this  mass  is 
already  small,  the  effects  of  compressibility  have  not  been  explored  any  farther. 
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APPENDIX  C 


VIRTUAL  INERTIA  OF  ENTRAINED  WATER  IN  BENDING  MODES 


Since  the  studie*tofXewis^  and  Taylor'^  on  the  effects  of  entrained  water  upon  ship 
hull  vibrations  made  practical  the  predictions  of  natural  frequencies  by  calculations/"'  ^  ^  'nc^ 
there  have  been  only  moderate  advances  in  the  understanding  of  water  inertia.  Because  the 
complexity  of  the  water  inertia  problem  it  has  been  difficult  to  include  these  advances  in 
the  calculation  methods  that  are  applied  to  ships. 

The  most  practical  method  of  including  the  water  inertia  in  the  calculation  of  ship 
bending  modes  still  appears  to  be  that  used  by  Lewis  and  by  Taylor.  This  consists  of  com¬ 
puting  the  water  inertia  per  unit  length  of  each  section  of  the  ship  as  though  each  were  an 
infinitely  long  prism  moving  in  parallel  motion.  The  result  of  these  calculations  is  a  curve 
of  sectional  virtual  masses.  These  sectional  inertias  are  then  reduced  by  a  factor,  generally 
designated  by  J,  which  represents  the  reduction  in  inertia  that  occurs  because  of  the  longi¬ 
tudinal  shape  of  the  hull  and  the  modal  pattern  of  vibration.  The  values  of  J  are  generally 
taken  to  be  the  same  as  those  that  can  be  computed  for  a  translating  and  a  vibrating  ellipsoid 
as  developed  by  Lamb'®.  Lamb'®  in  his  theoretical  preservation  includes  a  table  of  values 
of  J  for  a  translating  and  for  a  rotating  rigid  ellipsoid.  Both  Lewis  and  Taylor  determined 
J  for  the  bending  ellipsoid.  More  recently,  Macagno  and  Landwpber"  developed  J  values 
for  two  mode  and  three  mode  vibration.  The  shapes  of  the  modal  patterns  were  those  ob- 
tained  experimentally  on  the  GOPHER  MARINER^*  and  were  approximated  by  second  and 
fourth  order  equations  in  the  case  of  the  two  node  vibrations  and  by  third  and  fifth  order 
equations  in  the  case  of  the  three  node  vibrations.  The  J  values  that  are  obtained  for  pro¬ 
late  ellipsoids  of  different  ratios  of  major  to  minor  axes  that  are  vibrating  in  these  several 
modes  are  plotted  in  Figure  C-l.  The  surprising  thing  about  this  information  is  the  change 
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Figure  C-l  -  Longitudinal  Inertia  Coefficients  for  Motions  of  on  Ellipsoid 
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in  the  value  of  J  that  occurs  because  of  the  more  detailed  representation  of  the  same 
modal  pattern  as  the  approximating  function  is  changed  from  a  second  to  a  fourth  degree 
equation  in  the  two  node  case,  and  by  a  third  to  a  fifth  degree  equation  in  the  three  node 
case.  As  far  as  we  know,  there  have  been  no  determinations  of  J  for  higher  modes  of 
bending  vibration  than  the  3  node. 

The  effects  of  multi-noded  vibration  patterns  upon  the  J  values  are  bestjnvesti- 
gated  through  consideration  of  an  infinitely  long  circular  cylinder  vibrating  vertically 
with  a  sinusoidal  variation  of  amplitude  along  its  length.  This  problem  was  first  considered 
by  Kennard^  and  is  discussed  in  more  detail  in  Appendix  B  of  Reference  26.  According 
to  this  development 

i  _  k1  (a> 

k,  (a)  k^a)-  akQ(a) 


where  k  and  k,  are  the  first  and  second  modified  Bessel  functions  of  the  second  kind, 
o  1 

kj‘  is  the  first  derivative  of  k^  and  a  =  -jq  where  c  is  the  radius  of  the  cylinder  and  N  is 
the  distance  between  nodes.  A  plot  of  J  as  a  function  of  a  is  given  in  Figure  C-2.  It  is 
useful  to  have  this  function  expressed  as  a  simple  power  series  in  a.  Within  the  limits 
0. 6  <a<  4 


J  =  0.0609  +  0.7151 
a 


0.19093 


1 

1 

a 


and  within  the  limits  0<a<0.  6  J  1.030 


+  0.00297 

-  4  90 
a 


A  comparison  of  the  values  of  the  longitudinal  inertia  coefficient,  J,  as  determined  for  an 

ellipsoid  by  Macagno  and  Landweber  for  two  noded  and  three  noded  vibrations  and  as 

determined  from  Kennard's  relationship  is  given  in  Figure  C-3.  The  Macagno  and  Landweber 

study 

values  are  based  upon  a  detailed/of  the  potential  flow.  It  will  be  noted  that  the  infinite 
cylinder  simplification  matches  the  theoretical  value  of  J  quite  closely  for  the  2  node 
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Figure  02  -  Longitudinal  Inertia  Coefficients  for  a  Long  Cylinder  Deflecting 

in  Sinusoidal  Waves  of  Length  L 


CONESCO,  INC. ,  consultant!  In  engineering  science 


vibration  but  appears  to  be  low  (i.  e. ,  to  overcorrect  the  water  inertia)  for  the  three  node 
vibration.  Intuitively,  it  would  be  expected  that  the  infinite  cylinder  simplification  would 
match  the  higher  modes  more  accurately  than  the  lower  modes,  since  the  end  effects  are  of 
less  importance. 

For  the  present  hull  vibration  calculations,  it  was  decided  to  apply  a  longitudinal 
inertia  coefficient  correction  to  all  calculations.  The  magnitude  of  this  correction  is  taken 
to  be  0. 87  for  tow  frequencies  until  in  the  calculation  it  is  found  that  the  deflection  curve 
crosses  the  axis  at  two  points.  For  subsequent  calculations  the  longitudinal  inertia  coefficient 
for  the  calculations  at  any  calculation  frequency  is  taken  as  the  value  computed  for  the  next 
lower  frequency.  The  computer  ascertains  the  distance  between  the  most  remote  zero  crossings 
and  divides  this  distance  by  one  more  than  the  number  of  intermediate  crossings  to  find  the 
average  distance  between  nodes.  It  then  computes  the  value  of  a  and  the  value  of  J  for  a 
long  circular  cylinder  using  the  series  formula.  The  water  inertias  used  in  the  subsequent 
calculation  are  the  section  water  inertias  (2  dimensional  values)  multiplied  by  this  determined 
value  of  longitudinal  inertia  coefficient,  J. 
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APPENDIX  D 

STIFFNESSES  OF  THE  HULL  IN  THE  AXIAL  DIRECTION,  in  VERTICAL 
BENDING  AND  IN  VERTICAL  SHEAR 

The  sectional  stiffnesses  (stiffness  per  unit  length)  of  the  hull  are  computed  for  the 
axial  direction,  AE;  for  vertical  bending,  El;  and  for  vertical  shear,  KAG,  and  the  three 
values  are  plotted  in  Figure  2  along  the  length  of  the  submarine.  These  stiffness  plots 
reflect  the  discontinuities  that  occur  as  the  thickness  of  the  plating  changes. 

1.  Calculation  of  Axial  Stiffness 

The  main  contributor  to  the  axial  stiffness  is  the  plating  of  the  outer  and  inner 
shells.  Since  most  of  the  sections  are  circular,  the  cross-sectional  area  of  the  steel  is  easy 
to  compute.  This  area  is  corrected  for  the  effects  of  slope  and  for  the  stiffening  effect  of 
the  transverse  framing  (See  Appendix  E).  In  addition  to  the  shell  plating,  longitudinal 
framing,  keels  and  stringers,  where  continuous,  contribute  to  the  longitudinal  stiffness. 

The  contribution  of  such  heavy  longitudinal  structures  as  pressure  tanks  and  long  heavy 
foundations  to  the  longitudinal  stiffness  is  difficult  to  assess.  There  is  no  adequate  theory 
that  can  be  used  for  this  and  so  judgement  is  required.  In  Appendix  F  these  problems  are 
discussed  and  some  methods  developed  for  use  in  assessing  the  stiffness  contributions  of 
discontinuous  structures. 

Because  of  its  light  weight  and  its  flexible  connection  to  the  hull,  the  contribution 
of  the  superstructure  plating  to  longitudinal  stiffness  is  neglected.  Light  platform  decks  are 
not  included  for  the  same  reasons. 

The  stiffness  between  two  stations  is  found  by  forming  the  summation  of  l/AE  for 
the  sections  of  constant  area  and  finding  the  reciprocal  of  the  summation. 
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2.  Calculation  of  Stiffness  in  Vertical  Bending 

As  with  the  axial  stiffness,  the  major  resistance  to  bending  is  given  by  the  inner 
and  outer  shell  plating.  The  stiffening  effects  of  the  frames  is  included  in  evaluating  the 
bending  stiffness  but  not  any  correction  for  the  conical  shape  of  the  hull.  Discontinuous 
longitudinal  members  are  generally  of  less  importance  in  bending  than  in  axial  stiffness  and 
are  treated  in  the  same  manner.  The  center  of  gravity  of  the  section  inertia  is  taken  at  the 
centerline  of  the  ship.  This  is  not  entirely  proper  and  it  would  probably  be  advisable  in 
future  calculations  to  compute  the  location  of  the  center  of  the  structural  section  relative 
to  the  centerline  of  the  ship  and  to  correct  the  moment  of  inertia  for  the  offset.  Figure  D-2 
is  a  computing  form  developed  for  this  purpose. 

The  bending  stiffness  between  two  calculation  stations  is  the  reciprocal  of  the 
summation  (or  integral)  of  l/EI  for  sections  of  constant  moment  of  inertia  over  the  interval 
between  the  stations. 

3.  Calculation  of  Stiffness  in  Vertical  Shear 

The  stiffness  in  vertical  shear  for  a  unit  length  (KAG)  is  found  by  multiplying  the 
components  of  the  cross-sectional  area  by  shear  coefficients,  K,  and  the  shear  modulus,  G. 
Again  the  major  elements  in  the  shear  resistance  are  the  inner  and  outer  shells.  These 
elements,  since  they  are  circular  shells,  carry  a  shear  coefficient  of  1/2.  Other  elements 
that  provide  shear  stiffness  in  the  vertical  direction  are  vertical  sides  of  tanks  and  founda¬ 
tions.  For  these  the  shear  coefficient  is  one  and  is  applied  only  to  vertical  webs  working 
between  heavy  flanges.  |f  these  webs  are  discontinuous  the  transfer  of  the  loading  into 
and  out  of  them  must  be  estimated.  Horizontal  areas  do  not  contribute  to  the  shear  stiffness. 

When  the  overall  shear  stiffness  between  two  sections  is  required,  this  is  obtained 
from  the  sectional  stiffnesses  by  forming  the  sum  of  l/KAG  for  the  individual  sections  of 
constant  KAG.  The  reciprocal  of  this  summation  is  the  shear  stiffness. 
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CALCULATION  FOR  LOCAL  AXIAL,  SHEAR  AND  BENDING  STIFFNESS 
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Figure  D-2  -  Computing  Form  for  Determining  Sectional  Stiffness  Characteristics 
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APPENDIX  E 

INFLUENCE  OF  THE  CONICAL  SHAPE  AND  OF  THE  TRANSVERSE 
FRAMING  UPON  THE  STIFFNESS  OF  A  SUBMARINE  HULL 

A  conical  hull  has  a  somewhat  lower  axial  stiffness  per  unit  length  than  that  of 
a  cylindrical  hull  of  the  same  radius  and  hull  thickness.  A  hull  that  is  stiffened  by  frames, 
because  it  is  not  free  to  expand  and  contract  under  Poisson's  ratio  effects,  will  be  some¬ 
what  stiffer  than  the  unsupported  hull.  The  problem  of  the  conical  shell  is  considered  in 
Reference  27  and  the  effects  of  transverse  framing  upon  the  longitudinal  stiffness  is 
mentioned  in  Reference  28.  ( 

Consider  the  conical  shell  reinforced  by  circumferential  frames.  A  sketch  of  the 
shell  is  shown  in  Figure  El. 

V  i 

A 


V0 

Figure  E— 1  -  Diagram  of  Shell  Supported  by  Circumferential  Frames 
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Assume  that  the  effects  of  the  frames  are  uniformly  distributed*  along  the  length 
of  the  shell,  that  thfe  frames  can  carry  a  circumferential  load  but  no  axial  load  0(id  that 
the  ratio  of  the  frame  area  to  the  shell  area  per  unit  length  is  Let  the  pressure  normal 
to  the  axis  between  the  frames  and  the  shell  be  P^. 

The  following  equations  of  equilibrium  and  deflection  will  apply  (See  Ref.  27). 

r 

V  =  Vo  •  y-  El  (a) 

N  (the  meridional  force)  =  El(b) 

x  '  un<p 


Ng  (the  hoop  force)  =  rPr 


El(c) 


6  =  rt0 

E2(a) 

.x  r  X 

A  =  -  6  cot  <j>  I  +  -5 — r  I  e  dx 

sin  +  J  x 

v  1  X 

E2(b) 

A  ^ 

o  0 

c  u  ii  N0  Nx 

Forshell:  ee  "  eK  '  ►'EK 

E3(a) 

For  Frames: 

E3(b) 

Nx  N0 

cx  "  W  H  IT" 

E3(c) 

Equating  E3(a)  to  E3(b)  and  solving  for  N0 
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Therefore 


8  = 


+ 


1 

sin  <p 


( 


IllB 

T+T 


2 

) 


sin  <)>  sin 0 


V 

= - * 

Eh  sin  <p  (1  +  P)  J. 


(j  cos  ^+l+p(l-H  ) 


The  unit  length  stiffness  for  an  axial  load  is: 

_  2irr  h  Esin^  <£  (1  +  P) 

K - 5 - o 

JJCOS  (p  +  1  +  P  (1  “  (J^) 

The  correction  to  the  area  stiffness  is  therefore 

sin^>  (1  +  P) 

- T~ — -  2 

(i  COS  $  +  1  +P(1  “  p  ) 


For  a  circular  cylinder  <p 


IT 

1 


and 
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the  correction  factor  = 


1+  p(1-/) 


This  agrees  with  Schade's  value. 

The  maximum  correction  factor  will  occur  at  a  bulkhead  where  p  may  be  considered 
very  large  and  in  this  case  will  be  — =1.1  (about  a  10%  increase). 

1-F 

For  normal  framing  {3  is  about  0.4  and  the  correction  factor  in  the  cylindrical 
section  1.03,  i.e.,  the  stiffness  is  increased  about  3%  by  the  frame  support. 


For  an  unsupported  conical  shape  the  correction  factor  ^ — 

1  +  p  cos  <p 

For  the  steepest  conical  shapes  on  the  598  (the  ends  of  the  double  hull  section) 
<t>  =  66.  9°.  In  this  case  the  stiffness  factor  uncorrected  for  framing  =  0.  739  (about  a  25% 
decrease  in  stiffness).  For  the  outer  surface  of  the  hull  the  correction  is  normally  only  a 
few  per  cent  and  would  only  be  this  latge  near  the  bow  of  the  ship. 
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APPENDIX  F 

TREATMENT  OF  DISCONTINUOUS  STRUpTURES 


Any  complex  structure,  such  as  a  ship,  has  many  discontinuous  structural  elements 
in  it.  In  the  submarine  that  was  studied  the  more  important  discontinuous  elements  are  the 
tops  and  sides  of  pressure  proof  tanks  and  the  tanks  that  are  integral  with  the  foundations. 

The  problem  can  be  idealized  quite  easily,  since  the  hull  is  heavy  and  cylindrical 
bending  of  the  hull  can  be  neglected.  The  problem  becomes  one  of  the  stress  and  strain 
distribution  in  the  base  plate  and  the  stiffening  plate  where  neither  moves  normal  to  its 
plane.  The  three  examples  illustrated  in  Figure  F-l  cover  all  the  cases  necessary  for  the 
prediction  of  the  effect  of  the  discontinuous  structure  on  the  stiffness  of  the  hull. 

The  base  plate  can  be  defined  by  a  width  w  (=  2wr  for  the  cylinder)  and  thickness  h. 
The  discontinuous  plate  can  be  characterized  by  a  length,  I,  a  width,  b,  and  a  thickness,  t. 
There  appears  to  be  no  solution  to  this  problem  in  the  readily  available  literature 
and  although  it  would  be  possible  to  develop  a  solution,  this  would  require  some  time  and 
is  outside  of  the  scope  of  this  study.  The  results  of  similar  studies  of  the  effects  of  flange 
widths  upon  the  effectiveness  of  beams  are  presented  by  Schade  in  References  29  and  30. 

The  loading  condition  considered  by  Schade  that  would  most  closely  match  the  problems 
defined  above  is  the  beam  subjected  to  a  constant  moment.  However,  he  does  not  develop 
his  curves  for  this  case.  The  nearest  case  is  a  beam  carrying  a  uniform  load.  For  a 
practical  working  procedure  it  was  decided  to  utilize  the  curve  of  effective  breadth  ratio 
developed  by  Schade  for  a  flange  connected  to  double  webs  and  carrying  a  uniform  load. 

A  plot  of  this  curve  (Figure  10(a)  page  41  of  Reference  29)  is  given  in  Figure  F-2.  This 
curve,  when  entered  with  the  length  of  the  discontinuous  member,  L,  and  the  width,  B, 
equal  to  the  plate  width  when  fastened  on  two  sides  or  twice  the  plate  widthtyhen  fastened 
on  one  side,  gives  the  ratio  between  the  effective  plate  width  and  the  actual  pIcRewidth 
at  the  center  of  the  plate.  An  application  will  make  the  process 
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Fig.  Ft  Discontinuous  plate  attached  to  on  Axially  Stressed 
Plate  that  is  restrained  from  bending , 
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Figure  F-2  -  Effective  Breadth  Ratio  for  Plate  Attached  at  Both  Sides 
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clear.  One  of  the  ballast  tank  tops  in  the  GEORGE  WASHINGTON  is  20. 3  feet  long 
and  22.6  feet  wide  and  is  attached  to  the  hull  at  bottusides.  The  ratio  of  length  to 
breadth  is  0.  90.  Entering  the  curve  for  affective  breadth  ratio,  the  value  of  X/b  is  found 
to  be  0.  35.  Thus  the  effective  breadth  is  taken  as  0. 35  x  22. 6  x  12  or  95"  and  the 
stiffening  effect  of  the  tank  top  at  its  center  would  be  that  of  an  additional  steel  area  of 
the  hull  whose  area  is  95"  times  the  plate  thickness.  It  will  be  noted  that  this  stiffening 
effect  is  that  at  the  center.  As  the  ends  of  the  tank  are  approached  the  stiffening  effects 
will  fall  off  until  they  become  zero  at  the  ends.  The  rate  of  fall  off  is  assumed  to  be 
parabolic.  Thus  at  the  quarter  length  of  the  tank  top  the  assumed  equivalent  breadth  would 
be  3/4  of  95"  or  72". 
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APPENDIX  G 

PROCEDURES  USED  IN  DE-FINING  THE  SPRUNG  MASSES  IN  THE  SHIP 


One  of  the  primary  difficulties  in  determining  the  response  of  a  ship  is  the  definition 
of  the  structure  in  terms  of  mass  and  elastic  constants.  Normally  in  the  past,  the  ship  has 
been  defined  as  an  elastic  beam  to  which  masses  representing  the  weights  of  the  ship  structure, 
cargo,  and  entrained  water  are  considered  to  be  rigidly  attached.  More  recently,  attention 
has  been  given  to  representing  weights  that  are  not  rigidly  connected  to  the  main  hull  girder 
as  masses  connected  to  the  girder  tjhrough  a  simple  spring.  In  applying  this  procedure,  it  is 
necessary  to  define  both  the  main  hull  girder  and  what  constitutes  a  rigid  connection,  since 
obviously  all  connections  are  only  relatively  rigid.  For  the  purposes  of  the  vibration  calcu¬ 
lations  consider  that  vibrations  of  the  ship  are  those  of  the  main  shear  members  (usually  sides 
but  sometimes  also  bulkheads  connecting  the  upp4*  and  lower  flanges,  the  deck  and  the 
bottom)  and  that  the  rigidity  of  the  connection  of  all  weights  to  these  shear  members,  called 
the  main  hull  girder,  is  to  be  investigated  in  determining  whether  a  mass  is  sprung  or  not. 

Also  for  the  purposes  of  vibration  calculations  any  mass  whose  connection  through  the  ship 
structure  to  the  main  hull  girder  is  sufficiently  rigid  to  give  a  natural  frequency  if  held 
rigidly  at  the  hull  equal  to  1. 5  or  more  times  the  hull  frequency  that  is  being  investigated 
can  be  considered  as  rigidly  connected,  but  otherwise  should  be  considered  to  be  flexibly 
connected.  Flexibly  connected  members  can  be  concentrated  weights  (shock  mounted  equip¬ 
ment  are  of  this  category)  whole  portions  of  the  ship  structure  or  cargo  and  in  some  cases 
the  ship  bottom  with  its  water  inertia.  Because  of  the  importance  of  the  propeller  and 
shafting  system  as  a  flexibly  mounted  system,  particular  attention  has  been  paid  to  it  (See 
Appendix  H)  and  a  procedure  developed  for  finding  the  response  of  elastic  systems  (the 
hull  and  the  propulsion  system)  elastically  connected  (through  the  bearings). 

For  the  submarine,  fortunately,  most  of  the  structure, is  rigid  and  the  only  masses 
that  are  sprung  are  those  that  are  shock  mounted,  These  consist  of  definite  machine  components 
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generally  assembled  on  a  common  bed  plate  to  form  a  package  comprising  a  machine  assembly 
on  a  common  foundation  that  is  connected  to  the  hull  through  memberrthat  are  highly  elastic, 
the  shock  mounts.  The  characteristics  of  these  mounts  -  both  the  stiffness  in  several  directions 
and  the  damping  have  been  accurately  determined.  It  is  thus  possible  to  characterize  these 
sprung-masses  quite  simply  and  accurately.  In  order  to  set  a  reasonable  limit  upon  the  number 
of  masses  that  are  treated  only  those  masses  whose  weight  exceeds  2000  pounds  are  considered. 
The  items  considered,  their  weights,  stiffnesses  and  damping  and  location  are  given  in 
Table.  I. 

Because  of  the  large  mass  of  entrained  water  associated  with  them  and  because  they 
are  located  at  the  stern  of  the  ship,  where  the  effects  of  mass  are  pronounced,  the  most 
important  sprung  masses  are  the  control  surfaces  (stabilizers  and  stem  diving  planes)  in  vertical 
bending.  The  dynamic  responses  of  these  structures  to  a  harmonic  vertical  motion  was 
presented  in  Reference  2.  For  the  purposes  of  these  calculations,  these  surfaces  are  con¬ 
sidered  to  be  adequately  represented  by  the  sum  of  a  rigidly  attached  mass  and  a  simply 
sprung  mass  of  such  size  that  together  they  represent  the  predicted  response  of  these  structures 
which  was  presented  in  Reference  2.  In  future  studies  they  should  be  incorporated  in  the 
calculations  as  subassemblies. 
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APPENDIX  H 

PROCEDURES  USED  IN  DEFINING  THE  MASS- ELASTIC  CHARACTERISTICS 
OF  THE  PROPULSION  SYSTEM 


1.  In  Longitudinal  Vibration 

Although  the  propulsion  system  is  broken  by  a  dental  coupling,  the  longitudinal 
forces  that  can  be  carried  by  the  coupling  without  causing  slippage  are  generally,  larger  than 
the  longitudinal  vibratory  forces.  The  coupling  therefore  does,  not  isolate  the  propulsion 
system  from  the  propeller  and  shafting  in  longitudinal  vibration.  The  same  reasoning  applies 
to  the  dental  couplings  in  the  reduction  gear  and  the  turbines.  To  check  this,  consider  first 
the  Zur.n coupling  located  forward  of  the  emergency  propulsion  motor.  At  155  rpm  this 
coupling  transmits  241,  000  lb  ft  of  torque.  The  pitch  circle  of  the  coupling  teeth  is  23.  50" 
diameter  so  that  the  total  force  on  the  teeth  is  246,  000  pounds.  This  load  is  carried  by 
94  teeth  each  2-7/8"  wide  and  therefore  corresponds  to  a  load  of  910  pounds  per  inch  on 
each  tooth.  For  this  type  of  load  with  turbine  oil  lubrication,  a  minimum  coefficient  of 
friction  of  0.  10  is  reasonable.  Thus  the  coupling  could  carry  25,  000  pounds  before  Slipping 
would  occur.  This  is  well  in  excess  of  the  harmonic  force  across  the  coupling. 

With  the  couplings  considered  as  elastic  and  not  isolating  members  the  whole 
rotating  system  constitutes  the  vibrating  system.  This  comprises  the  propeller;  the  thrust 
bearing,  which  is  the  connection  to  the  hull;  the  propulsion  motor,  the  Zurp  coupling,  the 
reduction  gear,  the  turbines  and  all  the  connecting  shafting.  When  the  weight  and  longi¬ 
tudinal  stiffness  of  these  elements  are  studied  they  fall  naturally  into  a  number  of  discrete 
weights  connected  by  spring  elements  that  have  low  weights.  The  arrangement  is  shown  in 
Figure  5. 

The  weights  at  the  discrete  points  are  generally  formed  by  the  summation  of  manu¬ 
facturer's  weights  for  the  parts  but  in  some  cases,  for  example,  the  reduction  bull  gear,  it 
is  necessary  to  compute  the  weights  of  components  to  form  another  grouping.  The  weight 
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of  the  propeller  is  increased  by  a  water  inertia  for  axial  motion  with  no  rotation  as  given 
by  Lewis  and  Auslaender.  The  propeller  is  also  assumed  to  be  a  source  of  concentrated 
damping  whose  value  is  computed  by  the  methods  of  Reference  31. 

Proceeding  from  the  propeller,  half  of  the  weight  of  the  shafting  between  the 
propeller  and  the  thrust  bearing  is  assigned  to  the  propeller  and  ha|f  to  the  thrust  bearing 
mass.  Also  included  with  the  thrust  shaft  mass  is  the  weight  of  the  rotor  of  the  emergency 
propulsion  motor  which  is  closely  connected  to  it.  Because  of  the  longitudinal  flexibility 
of  the  web  of  the  low  speed  reduction  gear  (the  bull  gear)  the  next  natural  weight  division 
includes  the  bull  gear  shaft  and  hub  and  to  this  is  added  the  Zurn  disconnect  coupling. 

The  next  weight  item  consists  of  the  byil  gear  rim,  the  four  pinions  geared  to  it  and  the 
couplings  and  shafting  associated  with  the  pinions.  The  succeeding  weight  group  includes 
the  high  speed  gear  hubs  and  shafts  together  with  their  couplings  and  shafting.  The  rims 
of  the  high  speed  gears  with  their  mating  pinions  and  couplings  form  the  next  weight 
grouping  and  the  turbine  rotors  the  final  grouping.  The  weight  breakdown  of  these  items 
is  given  in  Table  HI. 

The  axial  stiffnesses  of  the  connections  between  the  weights  are  generally  computed 
from  the  lengths  and  diameters  of  the  shafting.  Where  the  shaft  is  pressed  into  the  propeller 
hub  and  into  couplings,  it  is  assumed  that  the  equivalent  elastic  length  of  the  junction  is 
one-third  of  a  shaft  diameter.  A  typical  shaft  stiffness  calculation  covering  the  section 
between  the  propeller  and  the  thrust  bearing  is  given  in  Tlable  H2.  A  very  flexible  element 
in  the  longitudinal  vibration  system  is  the  web  connecting  the  hub  with  the  rim  of  the  low 
speed  gear.  The  web  consists  of  three  plates  1 "  thick  that  connect  a  hub  that  is  about  31 11 
diameter  to  a  rim  that  has  an  inside  diameter  of  about  113".  The  web  plates  are  connected 
together  by  6  tubes  each  19"  in  diameter  and  on  a  72"  diameter  base  circle  that  have  a 
3/8"  wall  thickness.  The  webs  are  bored  out  with  18"  diameter  holes  in  many  of  these  tubes. 
They  also  incorporate  "sound  stoppers"  in  the  form  of  partially  flattened  circumferential 
tubes  of  about  4"  O.  D.  and  1/2"  wall.  Obviously  the  accurate  calculation  of  the  stiffness 
of  this  web  would  be  difficult.  As  a  guide  to  choosing  the  stiffness,  we  compute  the  stiffness 
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TABLE  H-l  -  Weight  Distribution  in  Propulsion  System  in 
Longitudinal  Vibration 


Station  A: 


Propeller 
Entrained  water 

1/2  Shafting  between  Propeller  and 
Thrust  Shaft 

Total 


23,550 

17,730 

10,400 

51,680 


Station  B: 

1/2  Shafting  between  Propeller  and  Thrust  Shaft 

Thrust  Shaft 

Emergency  Propulsion  Motor  and  Shaft 

10,400 
5,230 
16,  270 

Total 

31,900 

Station  C: 

Zuen  Coupling 

Second  Reduction  Gear  Hub  and  Shaft 

3,500 

10, 157 

Total 

13,657 

Station  D: 

Second  Reduction  Gear  Rim  and  Deadening  Rings 

Four  pinions  and  Associated  intermediate  Shafting 

11,010 

4,096 

Total 

15, 106 

Station  E: 

First  reduction  gear  hub,  coupling  and  shaft 

4  at  1000  pounds  each 

4,000 
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TABLE  H-l  (continued) 


Station  F: 

First  reduction  gear  rim,  web  and^ deadening  rings 

4  at  828  pounds  each 

3,312 

High  speed  pinions,  coupling  and  distance 

2  at  592  pounds  each 

1,184 

Total 

4,496 

Station  G: 

Turbine  rotor,  coupling  and  distance  piece 

2  at  3670  pounds  each 

7,340 
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TABLE  H-2  -  Typical  "Calculation  for  the  Axial  Stiffness 
"  between  the  Propeller  and  the  Thrust  Bearing 

Based  on  EB  drawing  9871-00142  Bu  Ships  No.  SSG(N)  598-203-1891801 


Main  Shafting  Details 


Section 

Length 

OD 

ID 

Area 

VA _ 

Prop,  to  Sleeve 

6" +6-1/4"* 

16.00 

4.5 

185.2 

0.066 

Sleeve  Section 

Sleeve  (bronze) 

75" 

16.00 

4.5 

185.2 

.  405  + 

.  396 

75" 

17.  373 

16.00 

9.5 

15.8  j 

Between  sleeves 

177-1/8" 

16. 00 

9.00 

137.5 

Stuffing  bx.  sleeve  sect. 

61" 

16. 065 

9.00 

138.1 

•441  ] 

Stuffing  bx.  sleeve  (bz) 

61" 

17.  373 

16.065 

8.9 

13.7* 

\  429 

Shaft 

91  13" 

91  TZ 

16. 003 

9.000 

137.6 

0.667 

Enlargement 

9-5/16" +  5.  7"* 

17.  25 

9.00 

170.1 

0.088 

Coupling 

16-1/2" 

28 

9 

552 

0.030 

Thrust  flange 

3" 

28 

616 

0.005 

Thrust  shaft 

23" 

15-3/4 

194.8 

0.118 

3.088 

*  1/3  of  diameter  for  fit 

+  f/A  multiplied  by  two  because  of  ratio  of  Elastic  Moduli  of  Bronze  and  Steel 


The  shaft  and  its  sleeve  are  in  parallel  so  (— )  , 

ct  comb 


1 

ct  shaft 


1 


ct  sleeve 

The  longitudinal  stiffness  between  the  propeller  and  the  thrust  collar  is 

AE_30xl06  0  _  in6  .... 

1 - OSH"  ~  9'7*  10  lb/,n- 

This  value  has  been  taken  as  10  x  10^  in  subsequent  calculations. 
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of  a  set  of  three  webs  loaded  by  a  unit  force  on  the  rim  and  restrained  at  the  center.  The 
deflection  would  be  that  of  a  circular  plate  rigidly  restrained  at  the  inner  and  outer  radii 
which  is  loaded  at  the  center  by  a  moment.  The  deflection  of  such  a  plate  is  given  by 

Timoshenko  as  <p  -  — .  Taking  an  inner  radius  of  15. 5"  and  an  outer  radius  of 

aaEh3 

56. 5"  the  ratio  becomes  0.  275  and  the  corresponding  value  of  a-  =  2. 7.  Since  the  moment 

p  ns  0  ns 

in  each  plate  for  a  force  F  is  x  — j-,  the  deflection  under  load  is  — x  </>  and  the  plate 
thickness  1",  the  deflection  for  a  unit  load  is  14.3  x  10"^  inches.  The  corresponding  stiffness 
=  0. 07  x  10  pounds  per  inch.  Because  of  the  stiffening  effects  of  the  pipes,  etc.,  this 
stiffness  was  arbitrarily  increased  to  0. 10  x  10^  pounds  per  inch.  It  is  not  anticipated  that 
this  stiffness  will  be  of  much  influence  on  the  hull  vibration  and  so  no  great  efforts  were 
expended  in  trying  to  rmprove  its  accuracy. 

A  very  important  stiffness  in  the  response  of  the  hull  to  the  propeller  forces  is  that 
of  the  connection  between  the  shaft  and  the  hull  -  1.  e.,  the  thrust  bearing.  For  this  reason 
the  stiffness  of  this  element  was  studied  in  considerable  detail.  The  calculation  for  this 
stiffness  presented  here  is  based  upon  summing  the  flexibilities  of  the  several  components. 

The  plans  that  are  used  are  the  following: 

Kingsbury  Machine  Works,  Inc.,  Philadelphia,  Pennsylvania 

Kingsbury  Dwg.  No.  464275  -  Bu  Ships  Dwg.  No.  F  3,  297, 494A 
37"  Kingsbury  Thrust  and  Journal  Bearing  Assembly 

Kingsbury  [>wg.  No.  464276  -  Bu  Ships  Dwg.  No.  F  3,  304,  231 A  -  37"  (6x6)  Kingsbury 
Thrust  Bearing,  Housing  Details 

Kingsbury  Dwg.  No.  464277  -  Bu  Ships  Dwg,  No.  F  3,  304,  232A  -  37"  Kingsbury 
Thrust  Bearing,  Bearing  Shell,  Oil  Catcher,  etc.  Details 

Kingsbury  Dwg.  No.  464278,  -  Bu  Ships  Dwg.  No.  F  3,303, 179A  -  37"  (6x6)  Kingsbury 
Thrust  Bearing,  Shoe,  Base  Ring,  Leveling  Plates,  etc.,  Details 

1.  Deflection  of  thrust  collar; 


The  bending  deflection  is  negligible 
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The  shear  deflection  is  given  by: 


c  _  3  P 

6s  -  s,  TO 


a 


1/2  loading  diameter 


13" 

7S 


•  r  =  1/2  root  diameter  =  1/2x15-3/4" 
a 

7  =  1.77 

h  =  plate  thickness  =  6. 375" 

G=  shear  modulus  =  11.8x10^ 

6s=  0.001812  x  10"6  in/lb 


2.  Deflection  of  Spherical  Surfaces  under  Shoes;  12"  spherical  radius  against  a  flat  surface 
6  shoes 

By  the  Hertz  formula  (Page  376  of  Ref.  33) 
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assume  F  =  120,000  1b  p  =  20,000  Rj  =  12"  E  =  30  x  106 


Then  ~ 

<t> P 


<PS 


=  0.8^/  ^  '  ]$•  -  - - 

(30)  x  10  x  12  x  20,000 

=  0. 137  x  10-6  in/lb 

=  0. 0228  x  10"6  in/lb 


3.  Deflection  of  upper  and  lower  leveling  plate 

Consider  this  deflection  to  be  the  sum  of  the  deflections  of  two  cantilevers.  The 

23" 

first  3.  8"  long.  4. 5"  wide  and  3"  deep  and  the  second  1 "  long.  4. 5"  wide  and  jg  deep. 
Since  there  are  six  leveling  plates  which  carry  the  load  on  two  sides  the  force  on  each 
cantilever  is  F/12. 

Therefore 


6 

F 


upper  leveling  plates  = 


1 

17 


p-sr  .. 

3  x  30  x  10*  x  4‘  5ty3 


x - t -  --3  =0.00  583x  10  in/lb 

3  x  30  x  10  x  4.  5  *(jj) 

U 

The  lower  leveling  plate  would  have  the  same  deflection. 


4.  Deflection  in  the  thrust  bearing  housing: 

Consider  as  direct  shear  between  the  center  of  the  shaft  and  the  flanges.  The 
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thickness  of  the  housing  is  1-3/4",  the  length  of  the  housing  carrying  the  shear  is  about 
30  inches  on  each  side  and  the  length  of  the  section  across  which  the  shear  is  transmitted 
is  about  20  inches. 

The  deflection,  8,  =  ^  yg  ]  ]  $  x  1 0^  x  10  ^  in/lb  .  Because  of  bending 

in  the  housing  and  the  flanges  this  is  probably  small.  Consider  the  housing  deflection  to  be 
0.03"  x  10'6. 


5.  Deflection  of  the  foundation: 

The  deflection  of  the  thrust  bearing  foundation  is  computed  through  the  use  of  the 
following  plans; 

EB  Div.  Dwg.  9841-057,  Bu  Ships  Dwg.  SSB(N)  598-1 1 2- 18906 16K  -  foundation  - 

Propulsion  Plant 

EB  Div.  Dwg.  9832-020,  Bu  Ships  Dwg.  SSB(N)  598- 101-1 890537F 

Frames  83,  84,  85,  86,  87,  88,  89,  90,  91  and  Tank  Flat 

EB  Div.  Dwg.  9833-006,  Bu  Ships  Dwg.  SSB(N)  598-1 14-189  0564C 

BHD  and  Framing  92,  93,  94 

The  thrust  bearing  sits  on  a  double  beam  about  5  frames,  i.e.,  130"  long.  Each 
beam  consists  of  upper  and  lower  flanges  about  20"  wide  and  3/4"  thick  separated  by  two 
webs  each  1/2"  thick.  The  spacing  between  the  flanges  is  about  22".  Therefore,  the 
moment  of  inertia  of  each  beam  is 

I  =  2x  15  x  II2  x  1  x<22>3  =4520  in4 
ii 


Assuming  that  these  beams  are  pinned  at  the  ends  and  are  subjected  to  an  end  moment 
ml 

the  rotation  will  be  Since  the  center  of  the  beam  is  about  36"  below  the  C  of  the  shaft 
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the  deflection  at  the  shaft  centerline  for  a  unit  load  will  be: 


6  F  x  36  x  130  x  36  -  ,ft-6  .  ... 

p  - - i  =  0.  208  x  10  m/lb 

r  3  x  4520  x  2  x  30x  10 

The  whole  foundation  structure  is  quite  loosely  tied  except  for  the  after  end 
attachment  to  the  tank  flat.  The  deflection  of  this  foundation  structure  in  the  longitudinal 
direction  under  a  unit  load  can  be  obtained  by  computing  the  Reflection  of  the  tank  top  to 
which  it  is  attached  in  shear  and  in  bending  and  adding  to  the  resulting  stiffness,  the 
additional  stiffnesses  of  the  longitudinal  members. 

The  deflection  of  the  tank  top  between  frames  89  and  92  in  shear  is  (since  there 
are  two  parallel  shear  elements) 


6  _  1  t 


1/2  x  50“ 


3/8"  x  78"  x  11.8  x  10 


=  0. 0725  x  10"6  in/lb 


The  deflection  of  the  tank  top  between  frames  89  and  92  in  bending  is  determined  as  follows; 

Because  of  shear  lag  the  equivalent  breadth  of  the  tank  closure  on  Frame  89  is 
obtained  from  Figure  F-2. 

X  =  165",  B  =  54"  ^  =  3.  g  =  0.86 

The  central  moment  of  inertia  determined  as  follows; 


Section 

Dim. 

Area 

C.  G. 

Mfr 

■fr  92 

Bhkd.  Fra.  92 

1-1/8"  x  120" 

135 

0 

0 

Tank  Top  web. 

1/2  x  78" 

39 

39 

1520 

19,800 

Tank  closure  Fr.  89 

1-1/8x54" 

61 

235 

78" 

20 

4760 

3255 

371,000 

39055 

169*000 

221,800 
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The  beam  is  loaded  as  follows; 


If  a  uniform  beam 

6  -5 x  VtJ  *  f2  *  5 
p  r 

=  - x 3 -  I  x  T  t  I 

30  x  10  x  222,000  L3U  J 

-  =0.0157  x  10'6 
P 

The  flexibility  of  the  tank  ^  >c> -  '32x10  in/lb 

The  shear  flexibility  of  the  longitudinal  members  are  as  follows; 


£ 

For  members  25"  off 


For  members  1 9"  off  C 

*  P 


_  _ 60 

60  x  ^  x  11. 8xl06 

=  0.  I7xl0'6 

-  60 

40X-JX11.8X106 

=  0.  25xl0"6 
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Since  the  tank  top,  the  two  longitudinal  members  25"  off  the  <£  of  the  ship  and  the  two 
longitudinal  members  19"  off  the  (  of  the  ship  are  in  parallel,  the  overall  flexibility  of 
the  foundation  connection  to  the  ship  is  the  reciprocal  of  the  sum  of  the  reciprocals  of  the 
parts. 


+  — »  +  - - 


aLtotal  0. 0882 x  IQ-6  0. 17xl0“6  0.25x10_6 


=  11.22xl06x  11.8x106  x  8x10_6 


o4oT  =0.0323x1  O'6  in/lb 


When  all  of  the  flexibilities  are  added  the  total  flexibility  of  the  connection  between  the 
thrust  shaft  and  the  hull  is  obtained. 


Element 
Thrust  Collar 

Spherical  Surfaces  under  shoes 
Upper  loading  plates 
Lower  leveling  plates 
Thrust  bearing  housing 
Rotation  of  found,  beams 
Axial  defl.  of  foundation 

Total 

The  corresponding  stiffness  between  the  thrust  shot 


Flexibility  in/lb 

0. 0018  x  10"J 
0. 0228  x  10“* 

0. 0058  x  10"* 

0. 0055  x  10“* 

0. 0300  x  10"* 

0.  2080  x  10"6 
0. 0323  x  10~6 

0.  3065  x  10"6 

and  the  hull  is  3.  28  x  10^,  say  3.  3  x  10^lb/in. 


2.  In  Vertical  Vibration 

A  sketch  of  the  parts  of  propulsion  system  that  participate  in  the  vertical  vibration 
as  a  separate  system  sprung  from  the  hull  is  shown  in  Figure  HI.  It  will  be  noted  that  this 


*9SH  r,  I 
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systam  consists  of  the  propeller  and  shafting  up  to  the  thrust  bearing  and  that  the  thrust 
bearing  and  the  elements  of  the  propulsion  system  forward  of  the  thrust  bearing  are  con¬ 
sidered  to  be  rigidly  connected  to  the  hull. 

The  weights  at  the  several  stations  on  the  shaft  are  those  of  the  propulsion  shafting 
between  half  way  to  the  station  forward  and  half  way  to  the  one  aft.  The  A^EI  is  formed  in 

the  same  manner.  The  propeller  weight  consists  of  23, 550  pounds  for  the  propeller  plus 

31 

4350  pounds  for  the  entrained  water.  Lewis  and  Auslaender  are  indefinite  as  to  the 

proper  value  to  use  for  entrained  mass  in  the  transverse  direction.  The  value  that  is  shown 

is  obtained  by  multiplying  the  value  for  entrained  mass  in  the  longitudinal  direction  with  the 

propeller  restrained  from  rotation,  an  entrained  mass  of  17, 730  pounds,  by  the  square  of  the 

ratio  of  the  transverse  projection  of  one  blade  at  0. 7R  to  the  axial  projection.  This  ratio 

=  pitch/0.  7wD  and  has  the  value  0. 495.  Thus  the  assumed  value  of  entrained  mass  in  the 

transverse  direction  is  4350  pounds  and  even  this  value  is  probably  high. 

The  moment  of  inertia  of  the  entrained  water  about  the  transverse  axis  is  formed  by 
31 

applying  Lewis  and  Auslaender's  suggestion  of  multiplying  the  value  for  longitudinal  water 

2  6 

inertia  by  D  /1 6  and  is  found  to  be  40.  9x  10  pounds  inches  squared.  The  moment  of  inertia 
of  the  propeller  itself  about  the  transverse  axis  is  taken  a  half  of  the  polar  moment  of  inertia 
and  is  therefore  18. 09x  10^  pounds  inches  squared.  The  sum  of  these  is  taken  as  59x  10^ 
pounds  inches  squared.  The  computer  program  does  not  treat  the  moment  of  inertia  as  con¬ 
centrated  with  the  mass  but  rather  spread  out  between  the  concentrated  masses.  In  accordance 
with  this  the  59x  10^  value  is  considered  to  be  distributed  between  the  stem  tube  bearing 
and  the  propeller.  A  more  precise  representation  of  this  inertia  would  be  obtained  by  treating 
the  propeller  as  two  masses  a  short  distance  apart  with  the  rotary  inertia  distributed  between 
them. 

The  damping  on  the  propeller  in  the  transverse  direction  was  obtained  from  Lewis 
31 

and  Auslaender  value  for  torsional  vibration  damping  by  the  following  procedures,  bg  by 
Lewis  and  Auslaender's  procedures  was  found  to  be  45,  300  ft  lb  sec.  It  is  assumed  that  the 
damping  torque  (=45, 300u8)  acts  at  0. 7  radius  (0.  7  x  *  5. 6  ft).  Thus  the  damping 
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45  300 

force  at  0.7  radius  =  -j —  can  be  related  to  the  linear  harmonic  velocity  at 

5,6  45  300 

0. 7R,  y  (=  5. 6u0)  to  give  the  net  damping  force  at  0. 7  radius  =  — - — j-  ^  =  1440  LTIb  ft. 

1 ' '  (5.6)  ’  sec. 

The  force  corresponding  to  a  transverse  harmonic  velocity  rather  than  a  radial  harmonic 
velocity  is  assumed  to  be  half  of  the  later,  or  720  If!  When  this  value  is  changed  to  inch 
units  and  multiplied  by  an  H  factor  of  0.60  the  assumed  value  of  vertical  damping  coefficient 
=  40  lb  in/sec  is  obtained. 

It  remains  now  to  determine  the  stiffness  of  the  coupling  between  the  shafting  system 
and  the  hull.  For  the  journal  bearing  that  is  incorporated  in  the  thrust  bearing  housing  and 
for  all  the  bearings  forward  of  this  since  they  are  oil  lubricated  bearings  and  are  mounted  on 
solid  structure,  the  weights  are  assumed  to  be  integral  with  the  hull.  The  elasticities  of  the 
connections  occur  in  the  stuffing  box  and  stern  tube  bearings  and  in  bending  of  the  shafting. 

The  flexibilities  of  the  bearings  arise  from  loco  ^deformation  of  the  hull  structure,  deformations 
of  the  bearing  support  structure  and  deformation  in  the  bearing  staves.  When  thefe  are 
evaluated  it  is  found  that  the  local  deformation  of  the  hull  is  negligible,  the  deformation  of 
the  bearing  support  structure,  a! though  designed  to  allow  some  accommodation  of  the  bearing, 
is  small  but  that  the  deformation  in  the  bearing  staves  is  significant.  The  stuffing  box  is 
assumed  to  give  no  support  to  the  shaft. 

The  flexibility  of  the  rubber  staves  in  the  stem  tube  bearing  is  determined  as  follows: 

The  shaft  diameter  in  the  bearing  is  17-3/8" 

The  length  of  the  stern  bearing  is  60" 

The  static  load  on  the  stem  tube  bearing,  taken  as  the  propeller  weight 
+  2/3  of  the  shaft  weight,  =  37,  400  lbs 

2 

The  projected  area  of  the  bearing  -  1042  in 

2 

The  load  per  unit  of  projected  area  =  35.  9  Ib/in 

From  Figure  H2  the  slope  of  the  deflection  curve  for  a  No.  10  stave  on  a  14-5/8" 

2  2 

diameter  shaft  at  a  load  of  35.9  Ib/in  is  0.001623  in/lb/in  .  Thus  the  overall  flexibility 
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of  the  bearing  is  0.  001623-f-  1042  =  1. 558 x  10  ^  in/lb.  When  the  flexibility  of  the  bearing 
support  tubes,  and  the  shaft  tubes  are  computed  they  are  found  to  be  0.  30x  10  ^  in/lb  and 
the  corresponding  stiffness  of  the  stern  tube  bearing  connection  between  shafting  and  hull 
is0.538xl06  Ib/in. 

For  the  forward,  stuffing  box,  bearing  the  projected  area  is  521  sq.  in.  and  the 

2 

load  is  assumed  to  be  7000  lb  giving  a  pressure  per  square  inch  of  projected  area  of  13. 5  Ib/in  . 

From  Figure  H2  the  slope  of  the  deflection  curve  for  a  No.  10  stave  with  a  14-5/8" 

2  2 

diameter  shaft  and  a  load  of  13.  5  Ib/in  is  0.00225  in/lb/in  .  Thus  the  overall  bearing 
flexibility  is  0.00225/521  =4.  31  x  10"6  in/lb. 

From  detailed  calculations  the  flexibility  of  the  forward  tube  and  the  bulkhead 
support  is  0.01  xIO  6  in/lb.  Thus  the  overall  flexibility  of  the  forward  bearing  is  4.  32x  10  ^1rj/lb. 
and  the  corresponding  stiffness  of  the  connection  between  the  shafting  system  and  the  hull  at 
this  point  is  0.  232  x  10  ^  Ib/in. 


H-17 


CONESCO,  INC.,  consultants  in  engineering  science 


APPENDIX  I 

GENERALIZED  BENDING  RESPONSE  CODE  (GBRC  1) 


The  generalized  bending  response  code  (GBRC  1)  can  be  used  to  calculate  the 
response  of  beam-spring  systems  to  specified  simple  harmonic  driving  forces  and/or  moments. 

For  example,  a  ship  hull  connected  elastically  to  other  elastic  systems,  say  a  propulsion 
system,  as  well  as  to  sprung  masses,  can  be  treated.  The  program  is  a  FORTRAN  II  program 
wh^ch  has  been  run  on  the  IBM  7090.  It  is  now  in  regular  use  at  David  Taylor  Model  Basin, 
and  is  referred  to  by  number  1-840-277-01. 

In  this  appendix,  the  difference  equations  used  in  GBRC  1  to  represent  a  non-uniform 
beam  (hull,  shafting,  etc.)  in  bending  are  derived.  Then,  these  equations  are  modified  to 
include  more  general  systems  such  as  systems  of  beams  connected  by  springs.  By  permitting 
complex  scaling  factors  for  various  input  parameters,  a  variety  of  damping  options  are  included, 
for  example,  Rayleigh  damping  and  hysteresis  damping.  Damping  coefficients  for  connections 
between  sections  can  be  given  explicitly. 

Longitudinal  or  torsional  vibration  problems  can  be  solved  using  the  program.  For 
example,  if  each  beam  is  made  up  of  sections  with  spring  connections  only,  the  equations 
reduce  to  those  for  a  mass-spring  system.  These  equations  are  those  used  to  approximate 
beam  systems  connected  by  springs  in  longitudinal  or  torsional  vibration  problems. 

The  program  contains  a  special  option  which  permits  the  effective  mass  for  each 
section  to  be  calculated  as  a  function  of  frequency. 

A  detailed  description  of  the  input  preparation  required  for  GBRC  1  is  included 
here.  In  Appendix  J  actual  input  sheets  for  the  particular  vertical  bending  calculations 
included  in  this  report  are  given. 

The  output  generated  by  GBRC  1  includes  a  tabulation  of  the  input  used  to  describe 
the  problem,  an  optional  summary  output  which  lists  amplitudes  for  deflections  and  moments 
versus  frequency  for  a  selected  set  of  sections,  and  finally  for  each  selected  frequency,  the 
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real  and  imaginary  parts  of  the  deflection  and  moments  as  well  as  their  amplitude  and  phase 
angle  for  each  section. 

Method  of  Analysis 

The  system  of  differential  equations  to  be  satisfied  by  the  hull,  shafting  system, 
etc.  in  bending  is  given  by 

(3. 1)  =  -  m(x)  -  c(x)  +  P(x,t) 


(3.2)  =  V(x,  t)  +  I  (x)ilbl  +Q(x,t) 

^  8t 


(O  2)  dY(*'  ^  -  ,r/y  t)  -  *) 

(3.3)  - r(x,  t) 


9y(x,  t)  _  M(x,t) 
9x  wr 


with  the  following  notation; 


x;  distance  in  the  longitudinal  direction  measured  from  the  origin  of 

coordinates 

t :  time  variable 

y:  displacement  normal  to  x  in  the  plane  of  bending 

y:  angular  displacement  relative  to  the  z-axis 

V:  shearing  force  in  the  direction  of  flexural  vibration  (y-direction) 

M;  bending  moment 

p:  effective  mass  per  unit  length 
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I  z:  effective  rotary  inertia  per  unit  length 

KAG:  shear  rigidity 

El:  bending  rigidity 

P:  external  forcing  function,  force  in  y-direction 

c:  damping  coefficient 

Q:  external  forcing  moment 

If  the  forcing  functions  are  assumed  to  be  a  simple  harmonic  function  of  time  so  that 

P(x,  t)  =  eiutP(x) 

Q(x,t)  =  eiwtQ(x) 

Then,  also, 

y(x,t)  =  e‘uty(x) 
y(x,t)  =  e'Uty(x) 

V(x,t)  =  eiutV(x) 

M(x,  t)  =  e'utM(x) 

and  equations  (1)  to  (4)  become 

(3.5)  =  H(x)  w2y{x)  -  ic(x)uy(x)  +  P(x) 

(3.6)  ^  =  V(x)  -  yx)u2y(x)  +  Q(x) 
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(3.7) 

(3.8) 


Solving  (3.6)  and  (3.7)  for  V(x)  and  y(x)  we  obtain 

2 

1  dM(x)  dy(x)  _(  lyz^W  \ /.a  QM 

+  ~3T  -V  fe(xT  +  KXSTxT 


or 


-I 


(3.9)  y(x)  =1 


/l  -  v*(x)  u 

KAG(x) 


1  dM(x)  ^ x  .  dy(x) 

KAGW  -xr  ■ QW  +  -$r' 


and 


or 


dM(x)  .  ,  i  \  2  dy(x)  w/  x  l  (x)  <■>  ^ 

”3x  +  VW“  =  V<x>  ^"tSgR  ) 


(3. 10) 


V(x)  = 


I  (x)w 
1  -  Jdill _ 
KAG(x) 


dM(x)  +  I  (xW2^ 

rar  +  3x 


We  now  subdivide  the  interval  from  x  =  0  to  x  =  L  by  the  abscissas 
xl/2  =  0<x3/2<x5/2<-"  <xN+l/2  =  L 


1-4 


CONESCO,  INC.,  consultants  in  engineering  science 


as  shown  below. 


1/2 


3/2 


E=i  \.  I  js . □ 

xN-l/2 


Then  if  (3. 5)  is  ingegrated  from  x  .,,tox  ...  .  ,  „  _ 

n-i/2  n+1/2,  for  n  =  1,  2,  3,  . . ,  N  ,  we  have 


V  i  *  -  V 
n+1/2  Vn-1/2 


I  n+1/2  x  , 

f  /  2  \  r n+ 

Mx)u  -  ic(x)ujydx  +  j  P(x)d; 

<ri— 1/2  Vl/2 


which  is  approximated  by 


(3.11) 


Vn+l/2 


V  wo 
n- 1/2 


2 

u-i(c  Ax)  u 
'n 


y  +  (PAx) 
n  n 


where 


(nA*)n 


(cAx)n 
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(PAx)  = 


■-J 


n+1/2 

P(x)dx 


n-1/2 


Now,  from  (3. 10)  we  can  approximate  for  n  =  1, 


N-l  by 


(I  Ax  ) 

f  n,  n+1 


(x)dx 
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and 


Xn  ~  ]//2  <Xn+1/2  +  Xn-l/2}  * 


The  end  conditions  to  be  imposed  determine  Vjy2  anc^  ^n+l/2  Part'cu'ar/ 


(3.13)  for  V(O)  =0  use  Vy2  =  0 

for  V(L )  =0  use  VN+]/2  =  0 


Substituting  (3. 12)  into  (3. 11)  and  taking  (3.  13)  into  account  gives  for  n  =  1,  2,  . N: 


(3.  1 4)  ot  n,  n+1  ^n+1  ^n)  +  ^n,  n+1  ^n+1 


-  at  .  (M  -  M  ,)"P  ,  (y  -  y  ,)  =  6  y  +  P  +  aC  ,  Q 
n-1,n  n  n-1  rn-l,n'n  'n-1  n  n  n  n,  n+1  n, 

■  ol  .  Q  . 
n-l,n  n-l,n 


where 


n+1 


['-( 


I  w 
.El... 
KAG 


-1 


0  otherwise 


'n,  n+1 


tt — r  for  n  —  1,  2,  3, 


N-1 


n,  n+1  ) 


[  V}  “  \n+l  f°r  n  =  !»  2/  3, . N-1 

'  r  n,  n+1 


0  otherwise 


n+1 
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2 

8  =  (u  Ax)  u  -  i(c  Ax)  u  for  n  =  1,  2,  N 

n  ^  n  n 

Now  (3. 8)  can  be  integrated  to  give  for  n  =  1,  2,  . . . ,  N 

f  Xn+l/2 

yn+l/2“yn-1/2  '  J  dx 

Vl/2 


which  can  be  approximated  by 


(3,15)  yn+l/2  "  yn-1/2 


where 


dx 

El 


From  (3.  9), 


yn+l/2 


can  be  approximated  for  n  =  1,  2, . 


N-l 


(3. 16) 


yn+1/2  " 


1 


■1  n 


M  . .  -  M 
n+1  n 


(KAG  Ax  ) 


n,  n+1 


yn+l  "  yn 


n,  n+l 
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while  y1/2and  yN+1/2 


are  determined  by  end  conditions. 


In  particular  since  M(O)  =  0,  and  V(O)  =  0,  we  will  use  the  end  condition: 


Similarly  for  the  other  end  we  will  use 


M 

n 


=  0 


Substituting  (3. 16)  into  (3.  15)  taking  account  of  end  conditions,  we  obtain  for  n  =  2,  3, . N- 


(3’ 17)  £n,  n+1  ^n+1  *  Mn}"  \  n+l  (yn+l  '  yn  > 


-t  (M  -  M  .)  -  ct  .  (y  -  y  ,)  =  5  M 
n-l,n'  n  n-1  n-1,  nwn  n-r  n  n 


where 


* n,  n+l  IkAGT  n+l 

n,  n+l 


for  n  =  1,  2, . . . ,  N-1 


S 

n 


for  n  =  2,  3,  . 


N-2 
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For  convenience  (3. 14)  and  (3. 17)  can  be  rewritten  as  follows: 


(3.18) 


'  Pn-1,  n  yn- 


,  -  ct  ,  M  , 
1  n-1,  n  n- 1 


+  (P, 


n,  n+1 


+  5  +  p 

n  rn- 


n-1,  n;  '  n 


+  (ot  ^,  +  ct  .  .  . . 

n,  n+1  n-l,n)M 


Pn,  n+1  yn+l 


°^,  n+1  ^n+ 


+  tt 


n,  n+1 


Q 


n,  n+1 


—  ac  |  Q  ■ 

n-l,n  n-l,n 


(3. 1 9)  -  <*  y  *  —  c  .  M  ,  +  (at  ,+at,  )y 
n-1,  n  n-1  n-l,n  n-1  n,n-1  n-l,n' 'n 


+  (t  .,  +  c  .  +r  )  M  -  cl  .  •  y  . ,  -  c  , ,  M  . ,  =  0 
n,  n+1  n-1,n  ’n'  n  n,  n+1  'n+1  n,  n+1  n+1 

For  n  =  1  and  n  =  N,  the  equations  will  incorporate  the  end  conditions  specified. 

In  summary,  equations  (3.  18)  and  (3. 19)  for  n  =  1,  2,  . . .,  N-1  can  be  rewritten  in 
matrix  notation 

(3.20)  a"z  =  T 
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where 


(3.  21) 


A  = 


A11  A12  0 


0  A_-  A 


\ 


\ 


\ 


'23 

0  \ 

\ 

'33 

A34 

\  \ 

\ 

\ 

V 

\ 

\ 

A 

\ 

\ 

\ 

\ 

'  0 


\ 


\  A  v  A 

_0  -  N,  N+l  AN,  N 


(3.  22) 


where 


(3.24) 


z  = 
n 


r. 

Z1 

r2 

P2 

1 

1 

1 

ii 

t“- 

1 

1 

1 

1 

1 

— ► 

ZN 

PN 

y 

P 

'n 

,, 

n 

P  = 

n 

M 

0 

n 
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with  P  =  P  -  ct  ,  Q  .  +  et  . .  Q 

n  n  n-l,n  n-l,n  n,  n+1  n,  n+1 


and  for  n  =  2,  3,  ,  N- 1 


5  4*  S  ,  4*  B  , 

n  rn-l,n  rn,  n+l 


Al  +  cC  i 
h,  n+1  n-l,n 


°n,  n+1  +  0Vi-l,n 
^n  +  en-l,n  +  ‘n,  n+l 


A 


n,  n+l 


A  i  - 

n,  n-1 

and  finally 


^n,  n+1 

-A 
n,  n+1 

-CL  J.1 
n,  n+1 

E  n,  n+1 

"Pn-l,n 

°h-1,  n 

-A  . 
n-1,  n 

-c  . 
n-1,  n 

61  +  P1,2 

*1,2 

*1,2 

?1  + 

ei,2:. 

6N  +  PN-1, 

N 

*n-i,n 

r  n  *  ‘N- 
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where  jj,  and  jj  ^  are  chosen  large  enough  so  that  Mj  and  will  vanish  for  practical 
purposes. 

The  above  formulation  can  be  generalized  to  apply  to  systems  of  beams  with  special 
connections.  In  this  case  each  beam  can  be  subdivided  and  the  sections  numbered  as  for 
example  in  Figure  It  1. Then  a  system  of  matrix  equations  of  the  form 

(3.23)  Al  =  ? 

can  set  up  where  in  general 


A  = 


Al, 

A12 

A13 

_> 

Z 

A21 

<N 

CM 

'< 

A23 

Z 

CM 

< 

AAA  - A 

N1  N2  mN3  NN 


Here  the  A  are  2x2  submatrices  of  A  and  z  and  P  are  defined  as  in  (3.  22).  In  general, 
n,  m 

A  will  contain  non- zero  elements  whenever  the  section  numbered  n  is  connected  to  the 
n,  m 

section  numbered  m.  In  particular,  if  n  and  m  are  adjacent  beam  sections, 

■P„ 


n.  m 


n,  m 


maL 
n,  m 


-d 
n.  m 


n,  m 
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If  sections  n  and  m  are  connected  by  a  spring 


A 

n, 


m 


+k  +  iuc 
n,  m  n,  m 

0 


0 

0 


where  k  is  the  spring  constant  for  the  connection  and  c  is  the  damping  coefficient 
n,  m  iti 

for  the  connection.  In  general 


n,  n 


m=l 

m^n 


A  +  D 
n,  m  n 


where 


In  addition,  if  section  n  is  connected  to  ground,  A  must  be  augmented  by 

n,  n 


-k  -  ic  u  0 

n  n 

0  0 


where  k^  is  the  spring  constant  for  the  ground  connection  and  c^  is  the  damping  coefficient. 
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When  the  above  formulation  is  applied  to  the  set  of  beams  with  special  connections 
illustrated  in  Figure  Inland  a  unit  force  is  applied  at  Station  6,  equation  (25)  takes  the 
form 

Ax=  P 
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The  Effective  Moss 

The  effective  mass  (u  Ax)  for  section  n,  as  it  is  introduced  in  equation  (3.  11) 

n 

and  later  used  in  the  definition  of  5n  in  equation  (3. 14),  is  made  up  of  two  parts.  The  first 
is  the  mass  of  the  ship  section.  The  second  is  a  frequency  dependent  virtual  mass  which 
takes  into  account  the  inertial  effect  of  the  surrounding  water.  The  method  used  to  calculate 
the  virtual  mass  is  that  discussed  in  Appendix  C  of  this  report  which  is  based  on  that  con¬ 
sidered  by  E.  K.  Kennard  (  26  ). 

This  code  provides  for  the  automatic  calculation  of  the  virtual  mass  coefficient  J 
as  a  function  of  frequency  for  each  section  n  of  a  system  and  then  the  effective  mass  for 
section  n  is  calculated  using 

(u  Ax)  =  m  +  J  (u)  m 
'n  n  '  '  n 

where 

m  is  the  specified  mart  of  section  n 
n  r 

mn  is  the  specified  water  inertia  associated  with  section  n 
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J(u)  is  the  virtual  mass  coefficient 

For  each  system,  the  displacement  of  each  section  is  calculated  for  a  set  of  frequencies. 
Until  the  frequency  is  such  that  the  real  part  of  the  displacement  for  a  section  has  at  least 
two  zeros,  an  initial  specified  value  of  J  will  be  used  for  that  system.  As  soon  as  the  dis¬ 
placement  has  at  least  two  zeros,  then  let  xq  represent  the  position  of  the  first  of  them, 
and.x^  the  position  of  the  last.  An  average  distance  between  the  zeros  of  the  real  part  of 
the  displacement  can  be  calculated  using 


d 


x 

a 


1 


where  n  is  the  number  of  zeros.  Then 


k,M 
J  *  t-  w  - 


where 

we 

*=  T 

c  is  the  system  radius,  and  kQ  and  kj  are  the  first  and  second  modified  Bessel  functions  of 
the  second  kind.  J  is  approximated  for  .6  <.  <x  <  4  by 


.0609  + 


.715 

cC 


19093  J 

r  + 


00297 

3~ 

cC. 


and  for  o  <  a  < .  6  by 

1.030  -  1.56- 
c 
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Input  Preparation 

The  input  data  required  for  the  GBRC-1  program  is  that  required  for  setting  up 
the  system  of  difference  equations  described  fcarlLer...  ..  In  particular,  the  following 
kinds  of  data  can  be  included  for  each  case; 

1)  The  total  number  of  beam  sections.  This  must  be  less  than  80. 

2)  The  range  of  frequency  values  and  the  frequency  interval  for  which  response 
curves  are  required. 

3)  The  system  to  which  each  beam  section  belongs,  i.e. ,  whether  to  the  hull, 
shaft,  or  some  other  system. 

4)  For  the  hull  anchor  any  other  system,  for  which  the  J-factor  is  required  for 
the  added  or  virtual  mass  calculation,  the  following  data  is  required: 

a)  An  initial  value  for  the  J-factor 

b)  The  radius 

c)  A  water  inertia  term  for  each  section 

5)  The  mass  of  each  section 

6)  The  integral  over  each  section  of  the  reciprocal  of  the  bending  rigidity, 
i.e.,  (Ax/EI)n. 

7)  The  integral  from  the  mid-point  of  each  section  to  the  mid-point  of  the  next 
of  the  reciprocal  of  the  shear  rigidity,  i.e.,  (Ax/KAG)n  n+j  or  (A^/KAG)n  m 

8)  The  integral  of  the  mass  polar  moment  of  inertia  from  the  mid-point  of  each 
section  to  the  mid-point  of  the  next,  i.  e. ,  (I^Ax)^  n+j  or  0^z^*)n  m 

9)  The  spring  constant  and  damping  coefficient  associated  with  each  spring 
connection. 

10)  For  each  set  of  parameters,  for  each  system  a  complex  scaling  factor  can  be 
supplied.  This  allows  damping  to  be  introduced  in  various  ways.  In  particular,  a  complex 
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scaling  factor  for  the  ship  mass  corresponds  to  the  introduction  of  a  Reyleigh  damping 
coefficient. 

For  each  section  and  connection  a  system  number  between  1  and  80  inclusive 
must  be  specified.  This  number  is  used  to  select  the  scaling  factors  to  be  applied  to  the 
parameters  for  that  section  or  connection.  It  is  also  used  in  order  that  the  program  can 
select  the  set  of  sections  and  connections  which  make  up  a  given  system  so  that  the  zeros 
of  the  displacement  of  that  system  can  be  determined  for  the  calculation  of  the  virtual 
mass  coefficient  J.  Section  numbers  for  sections  which  are  connected  to  each  other  should 
differ  by  less  than  15. 

The  input  data  cards  for  GBRC-1  are  laid  out  as  shown  in  Figures  H2  to.  1-6  . 

The  contents  of  these  cards  are  described  in  detail  below.  In  general,  when  a  series  of 
problems  is  being  run,  only  those  data  cards  which  contain  changes  from  the  previous  case 
are  required.  The  numbers  punched  in  columns  3  and  4  of  each  card  identify  the  kind  of 
information  that  card  contains.  An  input  deck  should  begin  with  a  title  card,  and  the 
set  of  data  for  each  problem  should  begin  with  a  data  control  card  (columns  3  and  4  contain 
90)  which  gives,  for  each  type  of  data  card  included,  the  number  of  such  cards.  The  data 
cards  for  the  problem  should  follow  in  the  same  sequence  as  they  are  given  on  the  data 
control  card.  The  last  card  of  the  input  deck  for  a  set  of  problems  should  contain  99  in 
columns  3  and  4  if  no  additional  problems  ore  to  be  run,  and  98  in  columns  3  and  4  if 
a  complete  data  set  for  another  unrelated  problem  follows. 

Detailed  formats  for  each  type  of  card  are  given  below. 

1)  Run  Title  Card 

Columns  Contents 

7  to  12  (62HAA 

13  to  72  These  contain  whatever  title  is  to  be 

associated  with  this  run. 
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Figure  1-4  -  Data  Input  Card  (3) 


PROBLEM  NO. 
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Figure  1-5  -  Data  Input  Card  (4) 
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Figure  1-6  -  Data  Input  Card  (5) 
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2)  Data  Control  Card 


Columns 

Contents 

3  to  4 

90  (Type  of  data) 

5  to  6 

Columns  5,  6  contain  number  of  data  cards  of 

7  to  8 

type  given  by  columns  7,  8  which  follow  for 
this  case. 

9  to  10 

Similar  to  previous  four  columns  for  next  data 

11  to  12 

cards  for  this  case. 

etc. 

Case  Title  Card 

Columns 

Contents 

3  to  4 

10 

7  to  12 

1S2HAA 

13  to  72 

Titfe  for  this  case 

Option  Control  Card 

3  to  4 

20 

9  to  12 

Opl :  Added  mass  option  selector 

0  if  added  mass  as  specified  in  input  is 
used,  1  if  J  water  inertia  as  specified  in 
input  is  to  be  used  (See  Section  4) 

13  to  16 

Op  2:  Selector. for  A-matrix  and  P  vector  setup 
option. 

0  for  present. 

17  to  20 

Op  3:  Selector  for  edit  routine  to  be  used  at 
each  frequency.  0  for  present. 

21  to  24 

Op  4:  Selector  for  final  edit  routine.  0  for  present. 
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Columns 

Contents 

25  to  28 

Op  5:  Selector  for  A-matrix  routine.  0  if 

A-matrix  is  not  printed;  I  if  A-matrix 
is  to  be  printed. 

Edit  Control  Card 

3  to  4 

£1 

9  to  12  ' 

Section  numbers  for  those  sections  for  which 

13  to  16 

displacements  and  moments  are  to  be  tabulated 

17  to  20 

*  vs  frequency. 

21  to  24 
25  to  28 

• 

General  Data  Card 

3  to  4 

30 

7  to  8 

Number  of  sections;  this  can  be  at  most  80 

9  to  16 

Starting  frequency  in  CPS 

17  to  24 

Upper  limit  for  the  frequency  in  CPS 

25  to  32 

Frequency  interval  to  be  used  in  CPS 

7)  System  Data  Cards  for  Added  Moss  Calculation 


3  to  4  21 

7  to  8  System  number.  The  hull  sections,  shafting 

sections,  etc. ,  can  be  distinguished  from  each 
other  by  assigning  a  system  number  to  each  set 
of  sections.  The  range  of  these  numbers  can  be 
from  1  to  80. 

9  to  16  Radius  associated  with  system  in  the  units  in  which 

the  unsealed  distances  between  adjacent  section 
centers'  arp  given. 

17  to  24  Initial  J-value  for  use  with  this  system. 
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8)  Section  ParametenCards 

Columns 


3  to  4 
7  to  8 
12 


15  to  16 
17  to  24 
25  to  32 
33  to  40 

41  to  48 
49  to  56 
57  to  64 
65  to  72 


Contents 

43  or  44 

Section  number 

End  condition  -  normally  0. 

1  for  V(o)  =  M(o)  =  0 

2  for  V(L)  =  M(L)  =0 


System  number 
Ship  mass 
Water  inertia 

<Ax/KAG>n,n*l 


V*1 


n,  n+1 


P  if  columns  3  and  4  contain  43,  Q  ,  if 
n  n,  n+1 

those  columns  contain  44. 


9)  Scaling  Factor  Cards  for  Section  Parameters 


Columns 


Contents 


3  to  4  41  if  this  card  contains  real  parts  of  scaling  factor. 

42  if  this  card  contains  the  imaginary  parts  of  the 
scaling  factors. 

15  to  16  System  number  for  which  scaling  factors  on  this 

card  are  to  be  applied. 

17  to  72  Scaling  factor  for  each  parameter  are  given  in 

those  columns  containing  that  parameter  value 
on  the  section  parameter  cards. 
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10)  Special  Connection  Parameter  Cord* 


3  to  4 

53 

7  to  8) 

The  section  numbers  for  the  connection  -  n  and  m 

11  to  12) 

if  m  =  0,  this  will  represent  a  connection  from 
section  n  to  ground. 

15  to  16 

System  number  with  which  connection  is  to  be 
associated 

17  to  24 

k  Spring  constant  for  connection 

ri/  m 

25  to  32 

C  Constant  damping  coefficient  for  connection 

n,  m  r 

33  to  40 

Cn  Coefficient  for  frequency  dependent 

damping  for  connection 

41  to  48 

(Ax)n,m 

49  to  56 

57  to  63 

Note  that  in  general  either  k  ,  C  and/or  C  or  (Ax)  ,  zrrt* 

n,  m  n,  m  ^  n,  m  '  'n,  m  KAG 


n,  m 


and  (|  Ax)^  ^  will  be  supplied  depending  on  whether  sections  n  and  m  are  connected  by  a 
spring  or  are  two  adjacent  sections  of  a  beam.  It  is  important  that  the  system  numbers  in 
these  two  cases  be  distinct,  if  J  factors  for  the  virtual  mass  calculation  are  to  be  calculated 
for  the  beam. 


1 1)  Scaling  factor  cards  for  parameters  in  special  connections 


Columns 


Contents 


3  to  4 


51:  This  card  contains  real  parts  of  scaling  factors 

52;  This  card  contains  imaginary  parts  of  scaling 
factors. 
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Columns 

Contents 

13  to  16 

System  number  with  which  these  scaling  factors 
are  to  be  associated 

17  to  56 

Scaling  factors  for  parameters  in  corresponding 
locations  on  special  connection  Parameter  Cards. 

1^|  End  of  PotE) 


3  to  4  99 
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